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PREFACE. 






-> The subject of Hydrodynamics em- 

braces many of the most difficult prob- 
lems in the range of physical research. 
Although, at all times attracting the 
-^ attention of the greatest minds, it is only 
'^^ within little over a century past that 
much real progress has been made in 
the solution of the many and compli- 
cated cases presented by the ordinary 
phenomena of Fluid Motion. The names 
of Euler, Lagrange and Laplace in the 
last century, and of Helmholtz, Stokes, 
Thomson, Rayleigh and Kirchhoff in 
this, stand out preeminently as those 
that have done the most to advance the 
theory to its present position. The 
object of the following article is to pre- 
sent in a short space the more important 
: points in the Mathematical Theory of 
'; Fluid Motion, as it has been developed 
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by these investigators. It is a want 
severely felt by any one making a study 
of this subject, that there exists no sepa- 
rate and complete treatise on Hydrody- 
namics. 

It is a fact, I think, greatly to be re- 
gretted, that the men who do the most 
for the real advancement of science so 
seldom present to the world the result of 
their labors and extensive knowledge, in 
any other form than an occasional me- 
moir in a scientific journal, or in a com- 
munication to a learned society. There 
are, however, notable exceptions to this 
general rule, as witness: Maxwell's trea- 
tise on Electricity and Magnetism; Ray- 
leigh on Sound, Cayley's EUiptic Func- 
tions, and a few others. If some one 
would present to the public a treatise on 
Hydrodynamics, of the scope of those 
mentioned on other subjects, he would 
certainly receive the gratitude of all 
physical students, and confer a great 
boon upon the scientific world. 

T. C. 
BAiiTiMoiiB, AprUf 1879. 



Elements of the Mathematical Theory 



OF 



FLUID MOTION. 



The following paper contains the 
mathematical investigation of some of 
the cases of the motion of incompressi- 
ble, frictionless fluids. The results ob- 
tained are to be considered as applying 
only to this class of fluids, unless the con- 
trary is expressly stated. The paper is 
intended to be introductory to a treatise 
which I hope before long to be able to 
pubHsh. 

In a subject so difficult as Hydrody- 
namics, there is but little chance for the 
discovery of hitherto unheard of proper- 
ties of the quantities dealt with, so, in 
what follows, the reader will not look for 
much that is absolutely new in the way 
of fact, although the arrangement of the 



work and in many cases the methods 
employed are my own. 

The references to the original sources 
from which information has been drawn, 
are given in every case, and I trust that 
these references, together with the mat- 
ter contained in this paper, will prove of 
value to any one interested in the most 
difficult but beautiful problem of fluid 
motion. 

Of late years, much has appeared in 
different places upon the subject of 
Hydrodynamics, but, so far as I am 
aware, there is no general work either in 
the English, French or German lan- 
guages. The aim of this paper and the 
treatise which will follow will be to 
combine in one work, all of importance 
that has been written upon the subject, 
and so enable the student to forego the 
immense amount of research necessary 
in order thoroughly to inform himself 
upon any one bynch of the subject. 

The short section which appears upon 
the theory of the Potential, is principally 
taken from Clausius's work upon that 



subject. The references to theoretical 
mechanics are, unless otherwise stated, 
to Thomson and Tait's Natural Philoso- 
phy. Kirchhoflfs Mathematische Physik, 
and Clifford's Elements of Dynamic, have 
also been consulted. 

§ I. 

GENERAL EQUATIONS OF FLUID MOTION. 

Let X, T, Z denote as usual the com 
ponent forces acting at the point {x, y, z) 
of the fluid reckoned per unit of its mass 
— ^then denoting by p the density of the 
fluid we haye for the forces acting upon 
the elementary mass pdxdydz the ex- 
pressions 

^pdxdydzy Y pdxdydz, Zpdxdydz; 

Now for the fluid pressure acting upon 
one face of the elementary parallelopiped, 
^^Ji ^y^z we have p 6ydz,p denoting the 
pressure on unit of area ; upon the oppo- 
site face it is, neglecting powers of 6x 
higher than the first, 



-<yy<y.(p+|tf.) 
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Consequently the resultant force due to 
fluid pressure acting in the direction of 
the axis x is, 

— dydz^dx. 

The equilibrium of this portion of the 
fluid therefore requires that 

dxdydz-^ —pX6xdf/dz= o 

with similar expressions for the other 
pairs of faces. We have thus for the 
equations of fluid equilibrium, 

These three equations can, of course, be 
replaced by the single equation of equili- 
brium. 

dp =p(K.dx + Ydy + Zdz) 

when dp denotes the variation of pressure 
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corresponding to the changes dx, dy, dfe, in 
the co-ordinates of the point at which the 
pressure is estimated. We see from this 
.equation that the expression, 

Xdx+Ydjz + Zdz, 

is either an exact differential or capable 
of being made so by a factor. If the 
forces X, T, Z belong to a conservative 
system, that is, a system possessing a 
potential, or for which the above ex- 
pression is an exact differential, we know 

that 

dK dY '. 

-= T- = &C. 

ay dx 

But when these conditions, are satisfied 
the quantity X-dx-^-Ydy + Zdz is an ex- 
act differential, or in the case of a system 
of conservative forces we have without 
the assistance of any integrating factor 

Xdx + Ydy -h Zdz=:dB> 

when E is the potential of the forces at 
the point {x,t/,z). It follows from this 
that dp^ — pd'R, or that /> is a function 
of R, then for all surfaces for which R is 
constant jo is also constani^ i.6. the^prpss- 
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ure is constant over all eqtd-potential 
surfaces. From these equations of 
equilibrium we can pass directly to the 
equations of motion, by means of D'Alem- 
bert's principle. Call w, v, w, the veloci- 
ties of a particle of the fluid whose 
co-ordinates at the time t are x,y,z, thus, 

dx dy dz 

u=.-=-. v=-^, 10=-— y 
dt' dt' dt ^ 

let w', v', ^'5 denote the accelerations to 
which these velocities give rise, then in 
our equations of equilibrium replacing 
X, Y, Z by, 

X-w', Y-v', Ti-w' 
we have for the equations of motion 

S=.(Z-»', 



wben we have ©f course, 
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,^du du dx du dy du dz 
^ ^'dt'^^'di'^dy'di^dz It 

Id d d d\ 
Hdt^''^^'d-y-^''Tzr 

We may for breyity replace this oper- 
ator by Y^, and we have thus for the 
equations of fluid motion the following: 
dp I J>u\ 



Tz-pv'-^tr 



Concerning the operator which we 

have denoted by =-, it is important to 

observe that it relates to a particular 
particle and not to a particular point in 
space ; the velocities w, v, w, are functions 
oi X, y, z and t, and denote the velocity 
which any particle has when it occupies 

the position denoted by aj, y, 2 and -=- 
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dt denotes the increase of velocity of a 

second particle over the one originally in 

this position, which arrives at this point 

after the lapse of time dt, while on the 

T)u 
other hand, ^T^^ denotes the change in 

velocity of the original particle during 
this time. When the motion is very 

small, the terms w-^, &c. may be neglect- 
ed, and we would have, 

Dt^di 

To our equations of motion it is neces- 
sary to add one more, expressing the 
continuity of the fluid. This equation 
simply expresses the fact, that during 
any natural motion there can be neither 
annihilation or generation of matter, or, 
referring to our problem, that the 
amount of fluid in any space at any time 
must be equal to the amount originally 
contained in that space, increased by the 
amount which has entered it during the 
time which has been allowed to pass. 
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diminished by the amount which has left 
it during that time. 

Let a, b, o, denote the co-ordinates of 
any particle of the fluid at an initial 
instant, on, f/, z, denote the values of 
these co-ordinates at the time t ; now in 
order completely to specify the motion it 
is necessary to express these latter quan- 
tities as functions of initial co-ordinates 
and the time. Suppose further, that da, 
6b, 6c, are the edges of a small parallelo- 
piped of the fluid, as these are assumed 
to be infinitesimal, the flgure will remain 
a parallelopiped during the motion. 

We have now for the co-ordinates of 
the extremities of the edges meeting in 
the point a, b c, 

a-^6a, b, c, a, b + 6b, c, a, b, c + 6c 

at the time t the co-ordinates of these 
points will be, 

dx^ dy ^ dz ^ 

^"^tJ^^ y^'fJ^^ ^+;7:;^^ 

aa aa aa 

• • • 

• • • 

From these we arrive, by a simple geo- 
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Sadbdc 



metrical process, at the volmne of the 
parallelopiped, which is then at the time 
t. 

dx dy dz 
dot da! da 
dx dy dz 
dV TV db 
dx dy dz 
dS d^ dc 

or representing the determinant by ^, 
jd Sadbdc; hence by our definition of 
continuity we must have 

j=i, 

or in general if p^ and p denote the initial 
and final densities of the fluid contained 
in this portion of space 

which simply expresses the fact that the 
density of the fluid contained in this por- 
tion of space must vary inversely as the 
volume of the space. 

This equation is known as the integral 
equation of continuity. The form of the 
equation most generally employed, how- 
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ever, is that which expressed the fact that 
the rate of diminution of density bears 
to the density at any instant the same 
ratio that the rate of increase of the vol- 
ume of an infinitely small portion of the 
fluid bears to the same ^infinitely small 
volume at the same instant. The sym- 
bohcal expression of this fact constitutes 
the differential eqication of continuity of 
the fluid. 

Let the flow towards the inside of an 
elementary parallelopiped of the fluid be 
considered as positive, then the flow to- 
wards the outside will be negative. 
Representing as before the edges of this 
elementary parallelopiped by dx, 6y, 6z 
we have for the flow through the face 
dydz in the direction of x and during 

the time dt 

pdydzudt 

through the opposite face the flow will 
be during the same time 

— 6y dzi pu 4- -^ . Sx\dt. 

These together give rise to an increase 
of mass 
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with similar expressions for the other 
pairs of faces respectively, perpendicular 
to the axis of y and z. Then the total 
increase of mass is 

J. js \ (d.pu , d,pv , d,pw\ 

but this increase of mass is also given by 

dxdydz-J-dt; 

equating these values and we have for 
the equation of continuity 

dp d,pu d,pv d,pw __ 
dt dx dy dz ~ 

or for incompressible fluids simply, 

du dv dw__ 
dx dy dz 

It may not be uninteresting to show 
how this differential equation of continu- 
ity can be derived from the integral equa- 
tion. We have, denoting the mass of 
this elementary portion of fluid by m, 
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m'=.p 



dx 


dx 


dx 


da' 


dh' 


dc 


dy 


dy 


dy 


da' 


db' 


dc 


dz 


dz 


dz 


6a' 


6b' 


6c 



6a 6b 6c 



Differentiating this with respect to t we 

have 

^Bp dA 



or 



_Dp p dJ 

the quantity --t- -^ will be found by 

easy reductions to be equal to 

du dv dw 
dx dy dz 

and the equation thus becomes 

Dp (du dv dw\ 

from which we obviously obtain the 
forms given above. Particular forms of 
this equation for special cases are often 
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quite simple — as, for example — suppose 
the motion of the fluid to be wholly- 
parallel to the plane xy ; in this case we 
have simply 

du dv __ 

dx d%f 

V du dv 

but ■^-= — 7- is the condition that the 
dx dy 

expression 

udy—vdx 

be an exact differential; calling ii dW 
we have 

_dW ^__d_W 

dy^ "~ dx 

The quantity W is called the stream 
function, and all motion takes place in 
the direction of the curves ¥^= const. 
If the motion be steady, the lines W 
= const, will form a system of tubes in 
the fluid, which may be called the tubes 
of flow. A much more general simplifi- 
cation of the equations of hydrodynamics 
exists however for certain classes of 
motion. It is a fact, the discovery of 
which is due to Lagrange that if at any 
time the expression 
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udx 4- vdi/ + todz 

is an exact differential, it will remain so 
throughout the motion; that is, if at 
any time we have 

du dv dw du du dv 



•) —; 7~5 



= 



dx dy dx dz dy dx^ 

4 

these quantities will remain so through- 
out the motion. Eepresenting these 
quantities by ^, 77, d, we may express 
this fact in another manner, viz. if at any 
time the motion of the fluid be irrota- 
tional, it will remain so during the entire 
motion. In particular, if the fluid origin- 
ally at rest be set in motion, by a system 
of conservative forces or pressures, there 
will be no motion of rotation throughout 
the entire motion. The following proof 
of this theorem is that given by Sir Wm. 
Thomson in his paper on "Vortex Motion," 
Edin. Trans. 1869. As this proof does 
not depend upon the quantity ^, we can 
give it in general representing by go the 

integral / — . We have then from our 

equations of motion 
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da)=zXdX'\-Ydt/ -\-Zdz 

/T>u , 
—l^r- dx-{- 






Now, according to hypotheBis, 

Xdx-\-Yd}/-^Zdz=€rR, 
and obvionsly, 



( 



Bdx T>du T>dz\ 



or smce 



jydx dDx 

dot = c?R — Y^ (w<]?tf; +vdf/-\- wdz) 

4- (wc?w + vdv + t<?^wj) 

or representing w' + v' + 1^' by V, 

D , 

YT-waaj •\-vdy-\- wdz) = d(R + JV" — o>). 

Integrating this along any arc (12) 
moving with the fluid we have 
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-(E+iV-;r); 

If the arc be a closed circuit the second 
member of this equation vanishes and we 
have 

=- / {udx 4- vdy + wdz) = o, 

or this may be expressed by saying that 
the line integral of the tangential com- 
ponent velocity around any closed curve 
of a moving fluid remains constant 
throughout all time. The line integral is 
called the circulation, and the proposi- 
tion may be stated. The circulation in 
any closed line moving with the fluid re- 
mains constant. In a state of rest the 
circulation is, of course, zero ; therefore, 
for the assumed case of motion generated 
by pressures or conservative forces we 
have that the circulation is always zero, 
so that udx -h vdy + wdz is an exact differ- 
ential. Bepresenting this quantity by d(p 

we have 

d(p da) dcp 

dx dy dz 
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from which we have for the differential 
equation of continuity 

d^q) d^cp d^(p __ 
or simply 

The quantity q) is appropriately called 
the velocity potential, and the velocity in 
any direction is expressed by the corre- 
sponding rate of change of (p. We may 
just here observe one fact concerning (p. 
If <^ is a minimum at any point, i,e, if it 
increases as we go away from that point 
there must evidently be a positive expan- 
sion of the fluid from this point in all 
directions. Similarly if 9> be a maximum 
at any point. Then th_e motion is in all 
directions towards this point and there is 
compression of the fluid. 

If there be neither expansion nor com- 
pression of the fluid within the region 
bounded by a closed surface, the great- 
est and least values of the velocity 
potential in that region must be on the 
surface; for since there is no expansion 
or contraction, there can be no maximum 
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or TniniTmiTTi value within this surface. 
If, therefore,' the velocity potential is 
constant over the surface, it must be 
constant throughout the enclosed region, 
since its greatest and least values are 
now equal. In particular, if it is zero 
over the surface, it must be zero through- 
out the enclosed region. When the 
velocity potential exists, the equation for 
determining the pressure can be put into 
a very simple form, viz. 

integrating 
but 



^f = 1+„-....„-=f^v 



so that 



to 



=/f=«-riHT-) 



or for our assumed case of incompressi- 
ble fluids 
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p \dt '/ 

Another form of the equations of fluid 
motion due to Lagrange is worthy of 
notice here, though the forms already- 
given, or Euler's equations, are those em- 
ployed in general in hydrodynamics. 
Since the quantities ic, y, 25, are functions 
of the initial coordinates of the point, we 
have 

dp dp dx dpdy dp dz 



= ±::77H- 



+ 



da dx da dy da dz da 



from these we have 



dp 



dx 



1 



dp dy dz 

da! doH da 

dp dy dz 

W dV db 
dp dy dz 

dc^ dc^ dc 



but we have also 



^x 

dt'' 



- 'J-=X — ^7n &C. 

p dx 



&c. 



Substituting the values of 



dp 
dx 



from 
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these last equations in those giving the 

values of -/^... and we have the Lagrang- 
da 

eian equations of fluid motion, viz. 

\df \da\de \dx 

+ -1 T-.— Z J- T- + - :^=o 
( dt )da pda 

with two similar ones containing b and c, 
respectively. The reader can see that 
from these equations we can readily pass 
to the forms given before. Where the 
forces X, Y, Z have a potential and there 
is also a velocity potential, these equa- 
tions become 

d'x dx d^y dy d^z dz^d(R—(i)) 
df da df da dt* da"" da 

d^x dx d^y dy d^z dz __ djB.—at 
d?db^ di^ib'^ dedh" dh 

d'x dx d^y dy d^z dz _ c?(B— <tf) 
dt~dc ~df do dfdc'" dc 

Differentiate the second of these equa- 
tions with respect to c, the third with 
respect to b, and subtract the latter re- 
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suit from the former. It will be sufficient 
to examine only two terms of the result ; 
we have then 

did^xdx) d { d^x dx ) ^dx d* dx 
del dt' db]'^dbi dedcl^dbdi* do 

dx d^ dx 
dc dt* db 

d {dx d dx dx d dx\__d \dxdu 



dt \ db dt dc dc dt db \ dt\ db dc 

dx du 



dx du \ 
"^dc db ) 



&c., 



the remaining terms will be obtained by 
advancing the letters. "We have then a 
quantity which differentiated for ^ is = 0. 
Performing similar operations on the re- 
maing pairs of equations we arrive readily 
at the following equations, where C,, C„ 
C, denote quantities which are independ- 
ent of the time. 

j dx du dxdu) j dy dv dy dv ) 
{dadb"^ dbda) "^{dadb db da) 

{ dz dw dz dw 



\dz dw dz «w) p 
[dadb dbda ) * 
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J dx du dxdu) j dj/ dv dy dv \ 
( dc da da dc) { dc da dadc) 

^dzdw dz d'w\__^ 
dcda dadc ) *' 

j dx du dxdu) j dy du dy du \ 
{ db dc dc db) ( db dc dc db ) 

\dzdw dzdw)^^ 
'^{dblc "Ic Jb)^^' 

Let now w^, v^, w^ represent the values 
of w, t;, w for ^ = o, thus at this time we 
have 

w=w„ v=v^, w=w^ 

x=za, y=b, z=c; 

Substituting these values in the above 
equations and they reduce inunediately 
to 

db " ~da "" » 

da dc * 
dv^ du), 

dc" db "^ * 

We have further, since u, v, w are func- 
tions of aj, y, 2, 
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du du dx du dy du dz ^ 
da ~~ dxda dy da dz da 

If now in our determinant A we de- 
note the separate minors by A^, B^, F^y 

^ _dy_ dz^ ^ Jy_ dz^^ 
^ db dc dc da ' 

substituting now the values of —= — &c. 

in the above equations C,, C,, C, and 
noting these last abreviations we have, 
since for incompressible fluids J=l, 

K Sdv dw 



\dy dw\ idw du) 

'\dz'~dy\'^'^^\dx'^dz\ 

. r J ^_^^ I __d^,_dw^ 
"^ '[dy dz]'' dc dh 

k \dv dw \ \dw du) 

'Idi'^ dy J "^^M^ "^) 

Ady dz)" da 

k \dv (Zt^ I j dw du \ 

'{Tz" d^]^^^\d^""dz\ 

■*■ »My dz\^ dh 



du^ 
dc 



dv^ 
da 
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Bepresenting the quantities on the 
right hand side, as we appropriately may, 
by 2 S^, 2 7/^, 2 Co, and solving the equa- 
tions for the quantities within the paren- 
thesis, we have 



^= 


=^. 


dx 
da 


+ V> 


db 


+ c. 


(2c 


V- 


=^. 


dy 
da 


+ V. 


dy 
db 


+ C, 


dy 
db 


:= 


=^. 


dz 
'da 


+ V. 


db 
dz 


+ c. 


dz 
dc 



IS the quantities ^„, ?;„, C^, which are 
the initial angular velocities of the particle 
of the fluid whose co-ordinates at < = o 
a, b, r,, are = o, we have that S, rj, C must 
also be = o, that is, . we arrive again at 
the theorem that if there be no original 
motion of rotation in the fluid there will 
be none at any future time. It will be 
of interest to obtain the equations which 
were used by Helmholtz in his great 
memoir on vortex motion. These are 
simply obtained from our equations of 
motion. The flrst of these equations 
written out in full is 
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du) •»- du du du du s ^ 
ax dt dx dy dz 

Now supposing the fluid initially at 
rest to be set in motion by conservative 
forces and pressures from the exterior, 
the analytic conditions for this are 

d^ dY ^. 
dt/ dx 

Therefore differentiating the first equa- 
tion with respect to y and the second 
with respect to x, and subtracting we 
eliminate co and the impressed forces, 
and have 

d dv d du duido du) 
"^dx dt dy dt dxldx dy) 

j d dv ^ ^*^ I 4. 

I dydx^^chi dx) 

from this we have obviously 



D __du ^ dv 
dz dz 



{ du dy 1^ 
{ dx dy\ 



and remembering that for incompressible 
fluids we have 

du , dv , dw 

— -I- — . -f zno 

dx . dy 'dz 
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this becomes 



DC __ du p dv dw - 
D^ ~" dz dz dz 

And similarly 

DS du i: , dv , dw ^ 
.=- = -.- ^+-^ ^+ --C, 
Dt ax ax ax 

The principle of the persistehce of ini- 
tially irrotational motion obviously fol- 
lows from these equations. 



§ 2. 

THE POTENTIAL. 

• 

It will perhaps be as well to make a 
few remarks here concerning the theory 
of the potential. It is not the purpose 
in these pages to go into that subject 
with any degree of fullness, but as there 
are a few leading principles which fre- 
quently recur a brief statement and 
derivation of them may be of assistance 
to some readers. 

We have, already observed one fact 
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concerning the velocity potential, viz., 
that if it is constant over a closed sur- 
face containing a certain dejGuiite region 
that it will be constant throughout this 
region, and in particular, if it be =o 
over the surface it will be =o through- 
out the contained region. When udx 4- 
vdy -j- wdz is an exact differential we have 
seen that by making it equal to d(p we 
can replace the quantity 

du do dw 
dx dy dz 

^ dx' ^ dy' ^ dz' 

Now let (T denote any closed surface, then 

if V be the outer normal to this surface 

we have -^ d<T for the rate of flow 
dv 

outwards through the element d(T in unit 

of time, then the total flow outwards in 

time dt is equal to 

'd(p^ 



ff 



dv 

where the integration extends over the 
whole surface. If the space enclosed by 
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<T be full both at the beginning and end 
of this time we, of course, have 



// 



dv 



This is the equation of continuity for 
the whole region. Applied to the element 
of volume dr or dx dy dz this gives us 

d^o) d*q) d^cp ., 

If we denote by r the distance between 
any two points ar, y, z, and a, ft, c, ^. e. 

r«=(a;-a)« + (y-ft)» + («-c)" 

and by m a constant we see that the 
equation -^'^=0 is satisfied by 

m ' 
^ r 

or - is a particular solution of the par- 
tial differential equation. 

The general solution is a homogeneous 
function of the quantities ic, y, g of the 
degree ^, where i is any positive integer. 
It is also well known that to every solu- 
tion of the degree i then corresponds 
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one of the degree— (t + 1) expressed by 
(pi 

The expression cpi is a Solid Spherical 
Harmonic of the degree L The expres- 
sion obtained by dividing q)i by r* which 
will be a function only of two quantities, 
viz., the angles 6 and ^, is a Spherical 
Surface Harmonic of the same degree. 

The quantity cp is now the Potential of 
the mass m upon the point (a, y, z). At 
infinity the Potential with its derivatives 
vanishes, but is finite and continuous 
throughout the space except at the points 
in which the masses are found, i. e., for 
a;=a, y=^, z=c. Let now in the space 
under consideration, which is supposed 
to be continuously filled with masses, m 
denote a mass placed at a given point A 
and let r denote the length of a line 
drawn from this point to any other B — 
then we know that the attraction of the 
mass m upon the point B is given by— 
dcp m 
dr ■" r"* 

Now suppose the line AB drawn to an 
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infinite distance, and farther that the 
space which contains the masses m is 
bounded by a closed surface. The line 
A3B will cut the surface an even number 
of times. Suppose now a sphere of radius 
unity to be described with center A, 
and then let the line AB describe a coni- 
cal surface cutting the element dm from 
the surface of the sphere, and ds^^ ds^, &c., 
from the given closed surface. Let € de • 
note the angle between AB and the outer 
normal to the surface; then where the 
line issueg. from the surface cos e will be 
positive and where it enters cos € will be 
negative. We have now 

, r'dw 
d8= 



cos € 



The normal force on the element ds is 
given by R cos. € ds; but R cos. € ds is 
equal to 

/dX. dY dZ\^ ^ ^ 

Kd^^d-y^dzr^y^' 

and therefore for the whole surface 
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Now since R= —j we have 

r 

R COS. € ds=±. mdu). 

As the point A is within the surface, the 
line AB first issues from the surface, 
giving a positive value of mdoj/ after that 
alternate positive and negative values of 
mdiu which destroy each other, so that 
we have simply 

2 R COS. € d8=mdat 
and 

jyB, COS. £ d8=m/ydw=4:7rm 

Now 

4:7tm=z4:7r///p dx dy dz 

d^ dY dZ 

' ' dx dy^ dz'~ 
or 

This equation holds throughout the en- 
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tire space which is filled with the masses 
m. 

If we represent an element of the 
space tinder consideration by dt and the 
density by p we may write 

from which 

•^ dc ^ dc r 
and finally 

dz •^ r \zl •^ dc d 

The second of these integrals is evi- 
dently a quantity of the same kind as 
(p'j the first is the potential of a mass 
which is spread out upon the surface 
fd8 and giving the surface density 



p COS. 



(")• 



Suppose now that <p' be the potential 
with reference to the point (x, y, z) of a 
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mass that is spread out upon a surface 
giving the surface density /i ; we have 

Assume the rectangular axes so that 2. is 
normal to the surface; assume on 25 a 
point infinitely near the surface, and sup- 
pose a circular cylinder with radius P' 
and having the axis of z for its axis of 
figure to cut the surface ; suppose P' in- 
definitely small but infinitely large 
with respect to the ordinate 2. Let cp^ 
denote the portion of cp belonging to 
that part of surface included in the cylin- 
der; the remaining portion (p'—q)^ will 
not become infinite or discontinuous by 
% becoming either =0 or passing through 
zero. 

By neglecting infinitesimals we have 

or (p' remains finite and continuous if 
the point under consideration passes 



39 . 

through the surface, i.e,, if \^z*=-±z. 
Further. 

or since, with respect to z, P' is infinitely 
great, 

dz ^ v's' 

or if « be positive 

if z be negative 

But (p'—(p^ is continuous and also 
— ^ , ^* ; consequently when js: changes 

from positive to negative passing through 

dcp' 
zero ■-^- changes suddenly by the amount 

—4^//. Now as we have taken z in the 
direction of the normal this fact can be 
expressed as follows: Call the inner 
normal v^ and the outer normal v^ then 
we have 
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which may be called the charcicteristic 
equation of ^' at the surface.* 

Suppose now t^^it we have a surface 
over which a mass is so distributed as to 
give rise to the surface density that we 
have denoted by ja giving then the poten- 
tial <p\ We will use, for the present, the 
symbol ^ to express a general function, 
which satisfies the equation, 

d'cp d?cp d*(p __ 

And we will take now U to represent the 
potential in the space under considera- 
tion, V for what we have denoted by (p\ 
the potential of a surface over which a 
mass is distributed, as mentioned above. 
Now we will introduce a new quantity 
W, which we proceed to define: At every 
point of the surface that we have spoken 
of conceive normals (positive) to be 
drawn, lay off on these infinitely small 
lengths, and through these points conceive 

* Maxwell, Elec. and Mag., p. 82, Vol. I. 
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another surface to pass, the elements of 
which correspond to the elements of the 
first surface and consider that on each 
element of this new surface a mass is dis- 
tributed equal to that on the correspond- 
ing element of the first surface but of 
opposite sign. Represent by h the 
negative product of the density of the 
mass on the element d6 of the first sur- 
face, by dG\ then if (a, ft, c) denote the 
element d6 and W is the potential of this 
element at the point (a, y, 2); we have 

d): 






„. d^ d^ dr 

SmCe, — = -J- ;r- 

dv drdv 



we have 

W=— / —5- cos, i^v) 

We can also express this in another 
way. Conceive a sphere of unit radius 
described about (a, y, 2,) ^s center; also 
conceive acone on d(S as base and having 
(aj, y, 2) as vertex, to cut the sphere, the 
area of the included portion being c?^9 
then 
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dv — 

the upper or lower sign to be taken 
according as cos. (rv) is positive or nega- 
tive. The cosine can only change its 

sign by {rv) passing through-^ ; evident 

ly this is the case only when the point 
can lie on a tangent to the surface ; then 
supposing that cos. {rr) does not change 
its sign, we have 

W=/:F ^(^^ 

when the upper or lower sign is to be 
taken according as cos. {rv) is positive or 
negative. If the above condition is not 
fulfilled, the surface may be divided into 
parts so that each part can satisfy the 
imposed conditions ; then will W be given 
as the sum of the corresponding expres- 
sions for each part. In order to examine 
whether or not discontinuity occurs in the 
value of W, by the point to which it refers 
approaching indefinitely near the surface 
— coinciding with it or passing through 
it — ^we will choose the axis so that z is 
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normal to the surface, and assume a 
point on z indefinitely near the surface; 
now by exactly the same process as 
that before employed, we see that for a 
negative z we have 

W, = — 2;rA; 

and for a positive z 

■W,= +27r^ 

Wj corresponding to the small portion 
cut out of the surface by a circular cylin- 
der of indefinitely small radius P, which 
is nevertheless infinitely large as regards 
2. Now since W, is independent of 25, W 
cannot become infinitely great by z be- 
coming infinitely small; and since Wj 
suddenly changes by 4k7rk, W does so 
likewise. From our equation for W, we 
have for W,, since it refers to an indefi- 
nitely small portion of the surface for 
which k is constant, 

since 

--5 i , dr dr dr 1 dr 

*'=A/i» +^ and-7-=;j— "3^ = j-f > 

^ ^ dz dr dz r dz 
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W =27rA:i —7--, 7^=^ \ 

from this follows 

or, since 2' is negligeable with respect to 

dz — ^'f'^p 
The second member being independent 

ot z, - - for z=o i& finite and continu- 
al? 

ous. Thus we see that the potential W 

for this double layer is finite always, but 

changes suddenly by the-amoimt 4;rA? on 

the point to which it refers passing 

through the surface in the direction of 

V ; the quantity ;i-iSj however, finite and 

continuous. With a perfectly arbitrary 

dW 

co-ordinate system the quantities --j-, 

dW dW .,, . , ^ ,. 

--7- J -,- ) will m general suffer discon- 
tinuity, since A; is in general not constant 
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over all the surface ; if A; be constant the 
diflferential co-efficients will suffer no dis- 
continuity at 'the surface. 

Suppose that we have two functions U 
and V of x, y, g, which with their derivar 
tions are single valued and continuous in 
the space under consideration, which is 
bounded by a closed surface. We have 
the identical equations 

dx dx dx^ "" dx^^ dx / 
dVdV ^d'Y d/rrdY\ 



dz^^dz' "dzK dzl 



dy dy dy^ "" dy\ dy 

d^dY 
dz 

Add and multiply by the element of the 

space city we have by changing in the 

second member a volume into a surface 
integral, 

(A) /'I— — ^ ^ ^ ^\dr 
•^ \dx dx dy dy dz dzl 

This equation expresses what is known 
as Green's Theorem. By an interchange 
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of U and V, which from the nature of the 
functionB can be effected, the first mem- 
ber of this equation will not change. We 
will have then 

= /*(?T(VJ*U-UJ*V) 

If U and V are velocity potentiale this 
gives, 

In the preceding equation (a) suppose 
U=V, then it gives, 

or 

2T= — y V-j — co8.(T^aj) + -J- COS. (vy) 

+ — cos.(r2) Vdff 

An expression for the energy in the form 
of a surface integral. 
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''The irrotational motion of incom- 
pressible fluid in a simply connected 
closed space -T is completely determined 
by the normal velocities over the surface 
2*. If 2" be a material envelope, it is evi- 
dent that an arbitrary normal velocity 
may be impressed upon its surface, which 
normal velocity must be shared by the 
fluid immediately in contact, provided 
that the whole volume inclosed remain 
unaltered. If the fluid be previously at 
rest, it can acquire no molecular rotation 
under the operation of the fluid pressures, 
which shows that it must be possible to 
determine a function cp^ such that 
A^<pz=o throughout the space inclosed 

by 2*, while over the surface —^ has a pre- 
scribed value limited only by the condi- 
tion 

"By Green's theorem if A^cp^o, 

/AS* 
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the integration on the right hand side 
extending over the surface i^, that on the 
left hand side over the volume. Now if 
cp and (p'\-Aq) be two functions, satis- 
fying Laplace's equation, and giving pre- 
scribed values of -~ , then the difference, 

dv 

/S(piBA function also satisfying Laplace's 

equation, and making -r— vanish over 

the surface of Z Under these circum- 
stances the surface integral in the pre- 
ceding equation vanishes, and we infer 

that at every point of 2, — ^', —~j-^ 

— ^— must be equal to zero. In other 

words, A (p must be constant and the two 
motions identical. As a particular case, 
there can be no motion of the irrotational 
kind within the volume 2", independently 
of a motion of the surface." 

The line described by a point in the 
fluid, moving always in the direction of 
the resultant velocity, is, as has been 
mentioned for a simple case, called a 
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stream line. If (p denote the velocity 
potential, we have obviously for the 
differential equations of a stream line, 

dx __dy __ dz 
d(p d<p d<p 
dx dy dz 

These lines evidently cut at right 
angles, the surface ^= const. 

If the normal velocity at every point 

of this surface is equal to zero, we must 

have 

dcp__ 

dv ~~ 

and this equation will represent a surface 
which cuts at right angles the surface 
^= const., and is consequently made up 
of stream lines; such a surface is appropri- 
ately termed a surface of flow. The 
following properties of such surfaces, 
though not directly bearing upon the 
matter in hand, will possibly be of inter- 
est: 

Suppose that we have for the equation 
of a family of surfaces 

fix, y, z)=q, 
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q being a variable parameter by giving 
constant values to which we obtain the 
equation of each member of the family. 

Then we have for the direction cosines of 
the normal in the direction in which q 
increases 

COS. (VX)=:V -^,QO^(vy)=iV-^ y 

ax ^ ay 

cosAvz) = y-rr 
dz 

if N be the component of the flow norm- 
al to the surface, w, v, w being the com- 
ponents parallel to sb, y and z respective- 
ly, we have 

\ ax ay dzl 

Hence, if N be zero, there will be no flow 
through the surface, which may then^be 
called a surface of flow; we have then for 
the equation of such a surface, 

dq ^ da , da 
ax ay dz 
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If there be another family of surfaces 
whose parameter is q\ and these are sur- 
faces of flow, then 

dq' do* da' 
dx dy dz 

If we have still a third family of surfaces 
whose parameter is q'^ that are surfaces 
of flow, then again. 



dq" . dq' 
dx dy 



+ V 



dz 



eliminating w, v and w between these, and 
we obtain 

dq dq dq 



= 



dx ' dy ' dz 
dq' dq' dq 
dx ^ dy ^ dz 
dq" dq^ dq" 
dx ^ dy ^ dz 

which is only satisfied by making q"=. 
some function of q and q'. Suppose 
that we have only the two first of these 
equations viz. : 

da dq dq 
dx dy dz 
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dq' da' dq' 
dx dy dz 

By eliminatiiig u, v, w in turn from these 
equations we arrive at the following where 
(P is an undetermined function of q and 

\dy dz dz dyl 
\dz dx dx dz I 

\dx dy dy dx I 

When one of the functions represented 
by q and q' is known, it is possible so to 
determine the other, that 4> shall be = 
unity. The flow in the direction of the 
normals to these surfaces being =o, this 
flow can only take place along the surface, 
and the intersection of the two surfaces 
will be a line of flow or stream line. A 
tube of flow or a stream fllament, is a 
tube whose bounding surfaces are made 
up of lines of flow. If the two paramet- 
ers q and q' have a series of values given 
to them, they will form a double system of 
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surfaces dividing space up into a number 
of tubes, each of which will be a tube of 
flow. We will go further into the con- 
sideration of stream lines in another 
place. 

§ 3. 

PLANE WAVES. 

Having thus briefly stated some of the 
more simple of the general properties of 
the equations of fluid motion we will now 
proceed to examine some of the problems 
which present themselves most naturally 
to the student. The first case that we 
shall take up is that of the motion of 
water in plane waves when the excur- 
sions of each particle are very small. 

When a body of water originally in a 
state of rest is endowed with a wave mo- 
tion each particle of the mass has a mo- 
tion of oscillation or, describes a closed 
curve in such a manner as to cause the 
particle of water, after a certain definite 
lapse of time, to resume its original 
position on the surface of the wave. By 
wave, is to be of course, understood sim- 
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ply the forma which the water assumes 
under the action of the disturbing force. 
Plane waves are those in which the mo- 
tion of every particle is parallel to a cer- 
tain fixed plane, and they may be gener- 
ated by bringing a soHd body, e,g,^ a 
cylinder, in contact with the surface of 
the water contained in a rectangular 
canal of uniform depth, and in such a 
manner that the line of contact shall be 
at right angles to the length of the canal. 
Plane waves will be generated the in- 
stant the contact takes place ; these will 
travel along the entire length of the 
canal; impinge on the ends and return; 
we will in our problem, however, limit 
ourselves at first to the case of a canal of 
indefinite length and then need not 
take account of the phenomena at the 
ends. Our mass of water being then 
supposed, contained in a canal as de- 
scribed, we will now proceed to the 
mathematical examination of the waves 
generated by such a disturbing force as 
has been mentioned. 

Assume the axis of z vertical and posi- 
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tive downwards, the axis of x parallel to 
the length of the canal, and that of y at 
right angles to its sides; the origin 
being on the surface of the water at rest. 
Let now ic^, y^^ z^ denote the initial vaiues 
of the co-ordinates of a particle and let 
u, V, to denote the displacements which 
the particle undergoes in the directions 
of jc, y, and z respectively. For plane 
waves advancing in the direction of ar, 
we have of course r=o, and the equation 
of continuity assumes the form 

the values of u and w being 

__dq) _<i(p 
"" d7^ "" dz' 

We have now to express ^ as a func- 
tion of X, 2, and ^, and, from the nature 
of the motion, periodic with respect to 
the last. 

The simplest way of expressing this 
periodicity, will be by introducing in <p a 
factor which shall be a trigonometric func- 
tion of the time. We may assume (p in 
the form 
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-f-<y«^, .sin. ( 

when <y is a constant, and the forms of f 
and F are to be determined, the latter 
is easily obtained by the following con- 
siderations. We know that cp must not 
change if we increase t by the time of 
oscillation, denoting this by r, and q} 
becomes 

-^<5z ^, ^ sin. (tt^. 

Now in order that (p may be a periodic 
function of ^, we must have 

F(«+T)=F(0 + 2;r 

which gives by expansion 

tF'(OH-j^F"(^)+ . . . = 2;r 
from which since r is constant, 

F"=:F'"=, &C. . • =0 

'« 2;r, 



F{0 =/^ -^'mt-^t 



and qi now assumes the form 

+ <r2^ ^ sin. %7t 
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Now for the determination of /; sub 
stituting this value of ^ in the differen- 
tial equation of continuity A^q^—o and 
it becomes 



dx^ 



+ (rV=o 



Integrating this and we have for f 
the equation 

/= A sin. (ra: + B cos. ax 

when A and B are the constants of inte- 
gration. This gives us now 

_ ± a^l . . . _ \sin. lit 

<P= r, (A sm.<ya;4-Bcos.(raj) — t 

This obviously may be written in this 
form 



r 



a^ 

€ 



I a,sm. (Tajsm. — t + o,cos. <ra;cos. — 1\ 

^(T2/, . 2;r^ , . 27r\ 

+ ^ (o,cos.(rarsin. — ^ + a,sm.<ya;cos. t\ 

-(P^l . . 2;r . ^ 27r \ 

+^ ( afiHL, (Tajsm. — «+yOjCos. (Toxsos. — t J 

-<''«/ /? ^ • 2;r^ , . 2;r \ 

^^ I /?jCOS. <ya;sm. — ^ + a^sm. (Tacos. — < j 
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By supposing the constants a, 5, ... . 
positive we can, by making the proper 
ones vanish and establishing certain rela- 
tions among the remaining ones, obtain 
an expression for q) which shall contain 
as a factor the sine of the sum or the 
cosine of the difference of the quantities 

In 

— t and (Taj. It is of course desirable to 

T 

introduce the quantity x into the trigono- 
metric factor as the form of q) is alike 
unaltered if we increase t by the time of 
oscillation or x by the wave length. Then 
making 

and for a simple advancing wave making? 
we have 

We have here before spoken of q> as 
the velocity function — there is a manifest 
appropriateness in this case in calling it 
the wave function — a name that we shall 
adopt for the present. 
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The quantity 

I <Tz — (Tz\ 

is called the amplitude of the wave, and 
is evidently the maximum value oi q); r 
is the periodic time, or period, after the 
lapse of which the values of tp recur; 
and (Tx determines the phase of the wave 
at the moment from which t is measured. 
It is evident that if we have any number 
of wave functions (p', q)" . . . which 
satisfy the differential equations 

A*(p'=o^ /l^(p"=o, &c. 

that this sum must also satisfy the 

equation 

A^S(p=io 

or any number of wave functions may be 
compounded into one resultant by simple 
addition. 

Before proceeding to the general 
problem we will examine the simple case 
of only one wave function. From the 
value given above for q) we have 
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to 






From these we see at once that the dis- 
placements u and w satisfy the equation 
of an ellipse whose semi axes are 

/ az -(yz\ 
/ az^az\ 

These values of ^, u and w can how- 
ever, be further simplified by finding 
what relation exists between the quanti- 
ties a^ and a,. To find this relation, we 
proceed to examine the forces which act 
on any particle of the fluid; these are 
well known to be of the form 



de ' ^' dV 

and the elementary equations of motion 
thus become 

1 c^ _ d^u 

p dx ^' da? 

1 dp _ d^u) 
"ply"" dt^^^ 
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from these we have 

+ J 9^^ f + const. ; 

but 

^f' _ __d_ d^cp /27rY dip 
'^'di^ "^ ~~ dx dt^~\ T / do^ 
cPw _ d d^<p_/27rY dcp 
— 'dt^''^ "" dz dt^ \T / dz' 
therefore, 

p=pU—] ^H-^gj+const. 

The constant is evidently p^ the initial 
pressure, the first part of the second 
number of this equation being the in- 
crease of pressure at the time t above 
what it was at the initial instant. 

As we suppose ourselves limited to the 
case of very small motions, we may re- 
gard u and w as quantities of the first 
order ; then it is obvious from the form 
of the expressions obtained for these 
quantities that (T is of the same order — 
the other factors being in general of finite 
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magnitude. Then any terms which con- 
tain <yw, or (Tw being quantities of the 
second order, may, with reference to 
those of the first order of magnitude, be 
discarded. 

These considerations now enable us to 
determine p wholly in terms of the 
initial co-ordinates. To do this it is only 
necessary to write x=x^ + it, z=z^-\-to, 
then from what has just been said we 
have at once 






Substituting these in the expression for 
^ and this again in the equation giving 
the pressure and we obtain for the 
latter 

(/27r\2/ o-z ^(rz\ 

sin. ^-^t + ffx^j H- gz^ + gw i +jt>. 

or as it may be written substituting for 
w its value 
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HC:n<./\.r-) 



+»4.V"-)(^(t'+'") 

Now for the determination of the con- 
stants, we observe first: that particles of 
the fluid originally on the bottom of the 
canal, must necessarily remain there 
during the motion ; second, for particles 
of the fluid whose z co-ordidates are equal 
to zero, i.e., for particles on the surface 
of the fluid at rest, p must be a constant 
=p^. Let h denote the depth of the 
canal and the first of these conditions is 
evidently reached by making w=o for 
those particles for which 2=h. This 
gives us then 
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^( ah -(Xh\ . /27r \ 



and consequently 



<yh -ah 



from which follows, 
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-ah 
a, e 



Again, make 2=0 and p=p^\ this gives 
obviously, 



from which 




ah -ah 


/2n'- . 




€ —e 



ah -ah 

£ +£ 

Now making for brevity 

-ah 
aa^e =— a 

And our expressions for u and w become. 



u 



^-% -£ ^jsm.(-^+(rajJ 



w 



By introducing the wave length now we 
can determine the constant a. Revert- 
ing to the expression previously given for 
^, and for convenience retain the two 
» constants a^ and a, — ^this was 



.-' 
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Let I denote the wave length ; we know 
that (p will remain unchanged by writing 
for X the quantity aj + /, as this simply has 
the effect of transferring the origin of 
co-ordinates from one end of the wave to 
the other. This substitution gives 

In order that (p remain unchanged, we 
must clearly have 

In 
(r=— . 

If we call io the velocity of each particle, 
we have also 

l=.rw^ or r= — 



Substituting those in the expression 
obtained above for p^ and it becomes 
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crh -ah 



l>=i>o+m+^]^^^y^ 



e +f 
-ah 



az^ -az. 












£ 



'sin. — (tt'^+sCo) 

r 



which reduces easily to 



27r 2;r 



/>=P. +i0^-2/>^«| 2^ 2^ [ s"i- 



B 4- f 

27r/ 



-5(.^+.) 



The same substitutions give us for cp the 
value 
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€ +€ (ott + X) 

and in like manner we can obtain for u 
and w the values 

( 2;r,, . 2;r,, J 2;r 
w = — a -J — - {n-z) — (hr-z) V COS. — 

e ■\-e {att-{-x)y 

j 2;r., ,2;r,, J . 2;r 
«^= — a-j — j-(A— «)— fA-2)>- sin. y 

The value of w is easily obtained. 

ah -ah 

It/ =1-77 ^^""-^h-^h^n 

£ + e 

27r^ 2;r, 

£ — £ 



2;r, 2;r, 



£ + £ 

from this is readily obtained 
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tt>= 



f 



In 



e - e 






V. 






^''h 
-t'' 



€ + € 

In the discussion of these values for 
^,'w, w a,ndp lies the whole theory of 
the motion of plane waves in a perfect 
fluid. We will now proceed to an exami- 
nation of these quantities. Denote by z' 
the vertical ordinate at the time t oi s> 
particle on the wave surface whose other 
co-ordinates are x and y. 

dz' --a ( "2;r,, 2;r,^ ,) . 2;r 

^=^= -^i rj (^^^) T (^-'^ r"^'— 

e — f {u)t-\-x) 

integrating 

, ap (2;r,, , 27r,, J 27r, ^ , . 

^ =2^1 t(^^) T^^-^) f "^^-T ("'■^^)- 

Differentiating this expression with re- 
spect to X and we obtain. 

dz' ai Z7t^, K 27r^, .) . 2;r 

€ — € (wt + x) 
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This yanishes for the values 

which three are the values of x and t for 
the points of maxima and minima of the 
longitudinal section of the wave. Differ- 
entiating again 

cTz' ap { 2;r,, , 2;r,^ J 2;r 

€ — € {<ot-\-x) 

I SI 

This vanishes for aJ=T> and -r-, with 

4 4 

T St 

t=j, and :j^^. Consequently there are 

points of contra-flexure at J and J of the 
wave length. It is obvious, from these 
considerations, what the curve is. 

It has before been remarked that u and 
w satisfy the equation of an ellipse ; this 
is now 
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u^ 



a'|~(A-^)'-f(A-^)[" 



e + f 



«• 



aN-— (A-2)— (Ah2;)J. 

The [plane of the ellipse is vertical and 
its longer axis is in the direction of the 
motion of the wave. Suppose now the 
particle under consideration to lie very 
near the surface of the wave — ^that is, z 
is very small as compared with h; then 

the terms containing — 7-(A-2j) may ob- 

ir * 
vipusly be discarded, and the only other 
terms which remain in the expressions 
for the semi-axes of the ellipse will de- 
pend on 

2zr,^ , 2;r^ 2;r 2;r' 

The equation of the ellipse thus becomes 
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u* V* 



47r 4;r 

The equation of a circle whose radius is 

As 

Of course the same result would be ob- 
tained by supposing the depth of the 
fluid infinite. Thus for particles near 
the surface of a body of water of finite 
depth — or for particles anywhere within 
the mass of a body of water of infinite 
depth — ^the motion is in a vertical circle 
whose radius is given above. Suppose 
again that the wave has an appreciable 
length — say l=h ; then for particles veiy 
near the surface the semi-axes become 
very nearly. 

2;r -2;r , 2;r -2;r 

or the path of the particle is nearly circu- 
lar — the ratio between these quantities 
being nearly 1.000,007. 
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The lengths of the axis continuously 
decrease as z increases. This is ob- 
vious in the case of the yertical axis giyen 

by 

e — e 

for as z becomes larger the exponents in 
this quantity become smaller, thus causing 
the first term in the brackets to diminish 
as the second increases, and consequently 
making the total value of the quantity 
diminish rapidly. Take now the hori- 
zontal axis denoted by 

e + « 
differentiating this with respect to z and 
we have 

For z<h the second of thest^ terms is 
always greater than the first and, conse- 
quently, -^ is negative; but the incre- 



J 
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ment dz is supposed positive, conse- 
quently d^ is negative, or the axis de- 
creases as z increases, Le,^ as we pass 
from the surface of the fluid. For z—h 
this axis becomes =2a, and the vertical 
axis vanishes. That is, for particles of 
water at the bottom of the canal there is 
only a motion of translation backwards 
and forwards in lines of length=2a. For 
particles at the surface of the fluid and 
for h^l the horizontal axis is nearly 
=535. 5a, or the ratio between the lengths 
of the horizontal axis at top and bottom 
of the fluid is nearly 267.7. 

Bef erring now to the values of u and 
w near the surface, we have 

w= — Afi COS. z!L («*< +05) 

V 

— r* 

2;r 
«j=A£ sin. — {(ot + x) 

Differentiating these for ^, squaring and 
adding the results, and we obtain the ex- 
pression 
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(t)'-(tr-f-r»)' 

The quantity on the left hand side of this 
equation gives the square of the veloc- 
ity of the fluid particle in its circular 
path; this, as we see from the second 
member of the equation, is independent 
of the time, and is directly proportional 
to the radius of the circle ; but the radius 

depends uppn the quantity ' l z for its 

value, and this increases as z decreases — 
therefore, for particles near the surface 
of shallow water, we have the velocity x 
varies inversely as their depth. From 
this it is evident that the water at the 
top of the wave moves most rapidly for- 
ward, while that at the bottom moves 
most rapidly backward. In the expres- 

sion for m discarding the terms ; we 

have for the velocity of translation of 
particles near the surface of water of 
finite depth, or anywhere within the 
mass of a body of water of infinite depth 
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Substituting for w its value of— and we 

T 

find for T the value 

9 

Calling t' the time of oscillation of a 
simple pendulum of length ^ we have 



^ 



■=)/lt'. 



7C 

From the above value of to we see that 
the velocity of transmission of the wave 
varies as the square root of the length. In 
all cases, indeed, the velocity is nearly as 
the square root of the length, for the 
factor 

€ " € 









is nearly equal to unity. In the case of 
very shallow water, the velocity dimin- 
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ishes considerably — as the quantity just 
written decreases rapidly with h — ^vanish- 
ing as is obvious for A=a 

So far we have confined ourselves to 
a single wave, that is, to a single value of 
(p satisfying the equation J*<p=o, But 
we have seen that if there are several 
values of tp each satisfying this equation, 
that collectively they satisfy the equation 

In the case when the wave lengths are 
the same but the phases different, we can 
easily find the result of adding together 
the waves given by the functions q), <p^t 

• • • <Pi. 
The value that we have tdready obtained 

for (p may be written 

-■ ^"-je + e ' Jsin.(rH+a;+a«) 

where it is to be understood that a^=o and 
is merely introduced for future symme- 
try. Any other function <pi which satis- 
fies the equation A^q)—o^ may be written 
under the above conditions 
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{(ot-hX + Oi ) 

And it is not difficult to see that a sum- 
mation of these functions will give us 

^ A ( (T(h^) -alhr^) ) . 

{(ot + XJfW) 

When 

A* = -2 Oi* + 22aj 2ak cos.cr («*— i— a* ) 

and 

^a^sin. ffa^ 

^^' ^= JSaiCOS. (rai 
j=o 

or inversely 

- _- ^a^sin. (T oi 

!F=itan. i=l 



^ -2'a^cos. (Toj 

If any of the quantities (To, &xi.,=7r or 
{2n'\-l)yr a change takes place in the 
summation. Suppose (Toi =7t then q^i 
becomes -(pi and is subtracted instead of 
added to the other functions; but if (Toi 

= (2/1 + 1) ;r we have since (T - - 
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ai={2n + l)^. 

That is if the difference of phase is an 
odd multiple of half the wave length, the 
corresponding wave function is to be 
subtracted instead of added to the others 
in finding the resultant of the system. 

Suppose now that we have two waves 
of the same length and amplitude, but of 
different phases and moving in opposite 
directions; the wave functions are ob- 
viously 

adding we have 

expanding the (arigonoinetric factor and 
reducing by aid of the relations 

COS. a=2 COS.* n"~l» ^iu. a=2 sin ^cos. ^, 



u 
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we readily find this expression to become 

(wj^+g) COS. a («— o). 

From this we obtain by differentiation 
the values of the displacements 

Dividing the first of these by the second 
we have 

aih-z) -(T{h-z) 

this ratio is independent of the time, 
consequently each particle moves in a 
straight line the inclination of which 
varies with x and z. Since also this ratio 



COS. 
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is the same at any given point, no matter 
what be the time, the w^ve is a standing 
wave or has no progressive motion. 

Hence if there exist in the liquid two 
waves having the same length and ampli- 
tude but moving in opposite directions 
the result is a single standing wave in 
which the particles move constantly in 
right lines whose inclinations to the axes 
vary with x and z, 

Eeverting now to our values of u and 
w, suppose that u=o, that is, that there 
be no horizontal motion, this gives us 

cos, — (wi-f)=o 

I SI (2n + l)/ 

or ctf^+a=7 or -. or ^ 7—^. 

4 4 4 

That is, there is no horizontal motion at 
the nodes of the wave. The greatest 
horizontal motion evidently corresponds 
to 

(ot-^XzziO,^ , or I 

or the greatest horizontal motion is at 
the highest point of the crest and the 
lowest point of the trough of the wave — 
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and evidently the motions at these points 
are in opposite dir^tions — ^which we have 
seen before from other considerations. 

In like manner by making to=o we find 
that at the top and bottom of the wave 
there is no vertical motion. Also, that 
the greatest vertical motion is at the 
nodes of the wave. 

Similar results are obtained by exam- 
ining the equations a. There is no 
horizontal motion at the points when 

X — H = T^^ -T-but there is a maximum of 
z 4 4 

vertical motion; also, there is no vertical 

motion at the points where a?—- ^^ = Oj -^> 

or ^ but there is the greatest horizontal 
motion. 

Airy has shown in his treatise on 
"Tides and Waves," that if the channel 
is of variable depth or vddth, that waves 
of the nature just described, that is, 
waves caused by the simple oscillation of 
the particles, could not exist by them- 
selves, but require for their existence the 
action of certain exterior forces into the 



82 

nature of which it is not necessary here 
to go. Without going into a mathemati- 
cal discussion of the reflection of wayes, 
I will merely state that after impinging 
upon a wall the particles of the wave 
move up and down the surface through a 
distance equal to twice their previous 
vertical displacement, and the same with 
particles at a distance of half a wave 
length from the wall; particles at a dis- 
tance from the wall of one quarter of a 
wave length, merely vibrate in a hori- 
zontal direction. When a series of waves 
enters shallow water the period remains 
the same, but the velocity and wave 
length diminish; the front of the wave 
becomes steeper than the back, and con- 
tinues to become more and more abrupt 
until the top of the wave curls over the 
front and the wave breaks in surf on the 
beach. 

CYLINDRICAL WAVES. 

If we throw a pebble into a body of 
water, or if we simply bring a solid body 
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in contact with the water at one point 
we know that a series of waves is gener- 
ated which are circular in form, concen- 
tric and having their center at the point 
where the disturbance takes place. The 
waves thus generated are called cylindri- 
cal waves, and the line passing through 
the center of these circles and normal to 
the surface of the fluid is called the wave- 
axis, and evidently is the geometrical 
axis of the concentric cylinders. 

In the case of such waves as this it is 
evidently not* admissible to assume the 
displacement in any direction as equal to 
zero, there will clearly be motion in the 
direction of all these axes. Our axis of 
z will be assumed as having the same 
direction as in the foregoing section, and 
the axes of X and Y will lie in the sur- 
face of the fltdd at rest. Our equation 
of continuity will have the general form 

the displacements being of course given 

by 
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dtp d<p d<p 

~rfa;' " dy* '~ dz' 

The same remarks that were previonslj 
made concerning the form of (p will hold 
here, the waves being supposed to ema- 
nate from the wave aads, so we can write 
for q) the equation 

when r as before denotes the periodic 
time. If the wave axis be taken as the 
axis of z we have, r denoting the distance 
from this axis to any point in a plane 
parallel to the plane of x y, 

and we may write (p in the form, 
± <>'« ^ / \ sill ( 2;r ) 

We must now, as before, determine the 
form of f. Substitute this value of ^ in 
the equation of continuity, and it is 
easily f otmd to reduce to the form 
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Transforming this by the subetitution 

8 

r=- 
<r 

we obtain a known form 

This is a particular case of the more 
general equation, 

of which a particular solution is the 
Bessel's function J» (s) given by 






+ 



2(2w4-2) ^ 2.4(2wH-2)(2n + 4) 



«• 



|;i + 6)"^'"*-f 



2.4.6(2^4-2) (2/i 4- 4) {% 

For our case n=o, and the function 
J„(«) is a particular solution, viz. : 

JoW=-j 1-2- + 2\4' ■" 27476"'"^ \ 

This is easily obtained directly, calling 
/*, the particular solution sought assiune 

/'^=a + a,* 4- a^8^ -|- 6,5* 4- a^ 4- 



86 

Substituting in the differential equar 
tion, and we have 

O = a,«-l + (a + 4a J + («, + 9a,)5 

+ (a, + 16a>'4- 

f yom which we have 

^ a « « « 



4' * 64 2«.4J 

when z is of course any positive integral. 
This gives us then for our particular so- 
lution 

( s* 8* 8* ) 

/.=J.(«)=«jl-2.-2r4.+2r5r6'+ -• f 

Designate now by ¥,(«) the other par- 
ticular solution of the differential equa- 
tion, and for brevity write simply T, and 
J, instead of Y,(«) and J^i^)- Let now 
I denote a function of 8, then it is well 
known that we can write 

Sabetitatmg this in the differential equa- 
tion and it becomes 

/I ^dl\d2 dV2_ 
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Dividing this by --=- and integrating, 
gives 

ny 

log. «+2 log. J,+log. -^= const. 

or assuming the consi^o, and passing 
to exponentials 

from 'which 

rfT _ 1 

da ■" «J/ 
and 

and, by substitution in the equation 
giving Y„ 

Now from the value of J^ it is clear that 
the expansion of -rzr-^ can only contain 
even powers of s and we can thus write 

multiplying by da and integrating gives 
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Y.= J. log. «+ J.j A ^+B ^ +0*^ + [ 
or as it may be written for brevity 

The quantity E^ is the product of two 
infinite series each of which contains only 
positive integral powers of 8 and, conse- 
quently, according to a principle in the 
theory of the Bessel's functions can be de- 
veloped in a series of these functions 
and, moreover, as all the powers in J, 
and the other factor of E, are even, only 
the even BesseFs functions will appear 
in the development thus 

Y,=:J, log. 8+aJ*'\-bP+cJ*'{- 

The co-efficients a, by c have to be determ- 
ined. 

Take again the differential equation 

fiiad perform the operation 

*'■'" s ds"^ 
on the quantity J, log. s and we find 
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icT Id ,)^ , 2 



2^ 

s ds' 



Hepresent the operator for brevity by J, 
then this is 

We have now from the general differen- 
tial equation affording Bessers functions 

and for n=o 
And for 7»>o 



AZn = Q^) *'^**~^ + J»+l f 



according to a known relation connecting 
these three consecutive functions. And 
80 we have, finally 

^ Jo log. 5= — J, 

Now from the above value of Y, we have 

^T,= J(J,log. 5) + a^J, 

+ ftJJ,+cJJ,+ 

And by aid of the transformations just 
;given 
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JY.= -^-J. + ^(J.+J.)+4'(J.+J.) 

+ y (J5+J,)x-- 

Now Y, being a particular solution of 
the differential equation ^/=o, we must 
have JY,=o/ this enables us to find 

b=-2c=Sd=-^=5(/=&G =2 

and by substitution 
Y,=J, log. 8+aJ^ 

+2[J-iJ.+4J-iJs+....] 

The complete solution of the differential 
equation Af^=^o^ i.€.y 

ds* 8 ds '' 
is now given by 

/=AJ,+BY, 

or 

=:(A+Blog.«) J.+BE, 

It maybe verified without much difficulty 
that the quantity E^ is of the form 

The quantities J, and Y, expressed in 
the form of definite integrals are — ^vide 
Boole's Diff. Equas. 



f 
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Jo=-y cos {s sin go) dco 

Y„=- / cos {s sin gj) log (4 8 cos'g?) (?c«7 
and we have for/ by substitution 

y= — / cos {s sin G?) (A+B log 

(4« cos* oi) d(tt 
or as this may be written 

/=/ cos (s sin tt>) (C + I>log {s cos'to)dfw 



when 

^^ A-hB log. 4 j^^B 

Before going on to the application of 
these results to the problem in hand, we 
will investigate the change produced in 
the quantities J^ and T^ by allowing s to 
become veiy great. 

Instead of / in the differential equa- 
tion ^f=o write //^ this equation thus 
becomes 

-^— ^r+4?)^^=^ 
and this for a very large is simply 
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— 5^— +/aA-o. 

This equation gives on integration 

f^=ia cos « + ft sin « 

or 

^_a cos 5 4-6 sin « 

when a and h are of course constants. 
We have then obviously from this 

_. __a cos «-h/? sin s 

^*- VT 

^af COS «+/?' sin 8 

from which we can see that for infinitely 
great values of 8 the functions J^ and Y, 
will vanish. 

Now by substituting for 8 its value of 
(fr we can, by taking as the argument of 
the functions thus obtained the quantify 

— 55, or J- write/in the form 

<=l(tl)' >=i; 



2;r 

+ cos 1 

T 
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For convenience of reference hereafter 
we shall write this in the form, 

/=AP(0)+B^(e). 

Substituting this value of /in the expres- 
sion for q) we pbtain, 

^^/'[AP (6) + Bn(0)] I sin ^^t 

\ 

or expanding this and writing instead of 
A and B, the quantities a„ b^, a„ 6„ a^, /J^, 

S' A» ^® ^^'^^ 

+Pcos-^(^^, +V / 
+/2sm— «(j,^^, j 

This gives us for w, v, w the following 
values ; 

2;r / <y2 -crgX) e?P 
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+ 









f . 2;r / Gz -(sz\ 
v= i sin — 1\„ e I of 

^ . 2;r^ / (Tz -€rz\ 

_ 2^^ / (TZ -(TZ\ 

. ^ . 2;r  / ixz -crz\ 

_ 2;r^ / ffz -<rz\ 

The expression for the fluid pressure 
is obtained here in the same manner as 
in the case of plane waves, and is 



w 



p=p\ (— ) <p+gz y + const. 
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Writing as before z=z^-\-w, remembering 
that <y and w are quantities of the first 
order of magnitude, and so discarding 
terms containing (Tw or higher orders, 
-we have 

P=P.+py^o'^P 1^ sm — « 
-t-Pcos — t 

{(t)'(A«"-)+H»?-v'')} 

+i2 sm — r 

{(V')'(I:^^Ar-)-Hv-.r')} 

^ 2;r 
+d2cos — « 

T 

For particles on the surface of the fluid 
at rest we have, of course, 2,=o and 
p=p^,' This gives us 

P sin. ^^t^(^) K4-a.)+M«,-«.)} 
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+/2cos^«{f-^)'|*.+A} 

In order that this may be satisfied we 
mnst have obviously 

«i+«i «8+«a ^+A *«+A ^<^ 

from which we obtain, 

^^^^^^ ^. 

We can now write 

. 2^- / az -cfz\ 

+ c,(r/2sin— «(^^^_^^^ ) 

^ ^TtJ +0-2 -(r2\ 



f 
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when the meaning of the constants 
c„ Cj, C3 is obTious. Introducing now the 
condition that w?=o for A=o we have as 

for plane waves 

-iSh 

a, ffh 

This gives us again 

(jh-crh 

. T / ^^ah -ah' 

B-\-B 

We have now for (p the equation, 
4p=a,Psin— «^^ -I- € / 

+ ^Pcos— <(/ + « / 

when a„ ^„ a„ ^„ are new constants, 
whose values are respectively, 

a, c^ 2^ p_^x 
Th' (TA ' ah' JJh' 

S B B B 

or those multiplied by any arbitrary con- 
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stant. This value of q) may be written 
in the form 

<p^/jM -<r{h-z)\ 

[(a,P + a,/2)8in — t + (*,P + i,/i)coB- t ] 
from which we obtain 

[(a,P + a,/2)sm ^t + (5,P + \n)coB^ t ] 
We have, however, 

d -dO'dx-^W' dxy^ 4: ) 

(T*x dF 



2 dff *''• 



consequently, 



ff*x/ff(A-e) -ff(h-z)\ U„dP, dDx 

(fy/«r(A-») -<r(A-«)\ ( / dP dn\ 

From these equations we see that the 
path of the particle is always ia a plane 
passing through itself and the axis of z. 
The expression for jo becomes now 

p^p.+Pffz^'\-2pg(T^ ya^jr + a^Djsm— t 



2n) e^€ 



2^ ) 



(TZ^ -<r2. 



(Xh -ah 

The values which have been obtained for 
the displacements and the fluid pressure 
afford the complete solution of the prob- 
lem under consideration. 

The results obtained are, however, very 
much modified in the cases where the 
particles are removed to great distances 
from the axis. We have already seen the 
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change produced in the fonction /in such 
a case, viz., this quantity becomes 
a cos ^ 4- ^ sin g 

or since (f is a constant 

^_ Acos (Tr+Bsin (Tr 

Vr 
We might have so transformed our first 

obtained value of/ that the infinite series 
therein contained should have proceed- 
ed according to ascending powers of -5 

and thus obtained the same result; this, 
however, would have been a difficult pro- 
cess. I 

The quantities P, and JQ are now given 
by the equations, 

_j sin or ^ cos ar 

Vr vr 

Substituting these in our value for q) and 
we have. 



•J (Oj sm ar+a^ cos<rr) sm —t 
+ (5, sin (Tr+b^ cos or) cos — t \ 
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for the simplest case of waves we can 
make a,=ft,=o and cti=^h=^a; thus 

'-^j?<n> .-'*"-"[<''+-) 

Differentiating this for r gives us the 
radial velocity of a particle, representing 
this by tf and we have 

rt^±.l<T(h-z) -(r{h-z)\ 

f /2;r, \ a . l^n ^ V\ 

•^ (TacosC — t-\-<5r\ +^ sm ( — < + (frl >• 

discarding term containing -, 

also 

The expression for je> also becomes, 
after making the proper substitutions 
and reductions, 

6 + € 
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Introducmg now the wave length I we 
have for great values of r 

1 J^ 

and also for a first approximation, 
111 






27f 

we have, as before, (T=-j-, and I = r oo. 



Thus (p becomes now, 
with as before, 



_ \lg €- € 

which for great depths or for particles 
near the surface becomes. 

The same deductions are to be made 
here as in the case of plane waves, riz., 
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that for particles any where within 
the mass of a fluid of infinite depth or 
near the surface of a mass of finite depth 
the velocity varies as the square root of 
the wave length. Write now a(y=^—a 
and collect all of our expressions: 

__ I a \ "lit., . 2;r,- . ) 










sm -t(^^ + ^I 



cos 



MO 






sm 



2;r 27r 



P=PM9^ir^P9-l= 



^ f — « 



Vz-o 2;r, 2;r 



£ + fi 



sm-j- ( w^ 4- rl 



w — 



/ 
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6 + € 

We see from these expressions that the 
amplitudes of cylindrical waves differ 
only from those of plane waves by the 

factor — -p=. — or in cylindrical waves the 

amplitudes varies inversely as Vr — and 
for particles very remote from the axis 
the amplitudes will vanish; whereas in 
the case of plane waves we saw that the 
amphtudes were always the same for the 
same depth. 

From the expressions for the displace- 
ments we have as before 



?rf(M?f(Mr 



+ 



w* 
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Now if A be finite we have as before for 
particles within the mass of the fluid — 
except near the surface— that they move 
in ellipses whose plane is vertical and 
passing through the wave axis, and 
whose transverse axis is in the direction 
of r. Also the axes of the elhpse de- 
crease as r^ increases and for particles in- 
finitely remote from the wave axis they 
vanish, or these particles are at rest. 
The axes also, as in the case of plane 
waves, continuously decrease as z in- 
creases, and for 2= A the transverse axis 

becomes — ■;=r- and the vertical axis van- 

ishes as it should do. 

If z be very small as compared with A, 
the equation of our ellipse becomes 

+ "Trr=i 



4:7t 4:7t 

J- z — — ^ 

A' e * A' £ ^ 

when 2;r, 

A=: =^6 

This is the equation of a circle whose 
radius is 
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2 



That is, for particles near the surface of 
a mass of fluid of finite depth, or for 
particles any where within the mass of a 
fluid of infinite depth, the motion is io a 
circle. It is shown as in the case of 
plane waves that this circular motion is 
uniform. In fact, all the results that we 
have obtained for plane waves are trans 
f erable into the corresponding results for 
cylindrical waves by merely multiplying 

by the factor — -=-• 

Suppose now that we have a series of 
n waves of the same wave length and 
amplitude but of different phases, start- 
ing from the same axis; let these waves 
be defined in the same manner as in the 
case of plane waves and we shall have for 
the resultant wave function 

when, sin(r(cu^ + r + ^) 

A* = 2; ttj* + 2 ^ aj 2 au cos crUi j^) 

i=o j=o =j-\-l * ' 
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and 






- __^ 2i aism (Toi 
!P^=- tan ;=i 

^Oj cos (To^ 

If the wave lengths are the same, but the 
amplitudes different by reason of differ- 
ent initial values of r^, the change in the 
form of these quantities is very slight ; 
they become 

^ la — d\ 

^°« "" [k-l k) 

and 

_^ 2 — = s"i o'fli 
!r= - tan <=i ^v/Vi 



<=~ 



JS— L cos <yaj 

Suppose now that we have two waves of 
the same wave lengths and amphtudes, 
but of different phases and going in 
opposite directions. The resultant wave 
function will be 
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2a { (r(h'-z)--(T(h-z)) 



sin 



which corresponds to a standing wave. 
Differentiating for r and z we have for 
the displacement of 7 and to, 

sin ^(^-^) 

sin a'(«>< + o ) cos Wr— ^ j 

The ratio — is independent of L and 
w 

we make the same deduction as before, 
that the particles move in right lines 
whose inclinations to the axes vary with 
z and r. 

Suppose that we have a series of par- 
allel wave axes, and that waves proceed 
from them having the same length and 
equal amplitudes but different phases. 
Let the wave functions be given as 
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^0 (ajhr-z) ^(r{h-z) 



) 



"^■'-iv.. 



sin (r(o)t + r(®) + a J 
^(^(h-z) -(r{h-z)\ 

sin (r(<ii^+H^^ + a, 
^ _^rn_ J^(A_^) -^{h^)l 

sin (T{Gi)t + r^^^ + «n ) 
adding these we have 

sin a{Gat + r" + aJ sin ff{cDt + rO) + a^) 



{ 



si n(y(&j< + r('*)+c^>) ) 
from which 



• • • • I 



i=n 
W 



i=o \€ - € J 

sin(r((w«+rW4-ai) 



{^ 



• • • • I 



sin (T ((»« 4- r(**) + fln ) 1 

The waves may evid^itly so move that 
at certain points the vertical displace- 
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ments shall be equal to jsero. We can 
determine these points by placing the 
trigonometric factor of to equal to zero; 
thus 

sin <T {wt+f<^)-{-a^) sin(y(<tf^+r<^) + a^) \ 

sin(r(w«H-r<»*)4-an) 
+ 7-= '=0 

This expression can be divided into two 
parts, one of which shall have for a 
factor sin ffaot, and the other cos aaot 

( cos (T (H<^> 4- qp) 
cos (T (r(i)+aj) 

c os a (rW + g^) ) 



sm <T GO t 





"T • • • • 



_ ( sin (T (/•(«)+ aj 

sin a (r(^) 4- «i 

I' • • • • 



]-' 



Equate separately to zero the factors 
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multiplying sin ff m t and cob (T to t, 
square and add tke resulting equations 
and we have after some easy reductions, 

2i — 77x+ ^ ^ COS 0< 7= f =0 

For the simple case of w=2 or two wave 
axes we have since a^z=zO 

—7-,+ --771 + —-== cos(r(r<<')-r(i)-^)=o 

If r,(*)=r/^) this becomes 

2 2 

— 4-— cos a (— a, )=0 

' ' 

If now -<Ta= + {2n + l)7t this equation 
will be satisfied, i.e., if 

Therefore if the difference of phase is an 
odd multiple of half the wave length the 
vertical displacement is zero — ^but only 
for the points for which r,<®)=r/^>. The 
points defined by the equation 

' 'o 

lie on a plane which from its relation to 
the waves may be called the plane of 
symmetry. We will now examine a little 
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more closely the conditions at this plana 
of symmetry. We have 

sin 0(oft+r) 
,_ a, j a{hr-z) -a{h-z) ) 

sin (y(aj^4-r' + aj) 

writing, for convenience, r' for r^^\ We 

have also 

_ dtp dr dcp' dr' 

~ dr dx dr dx 

jA(pdr d(p' dr^ 
~drdy dx' dy 

when r^^^-^y" and r'*= (a— 2a)" + y'. 
Now since 

we will have for r—r* 

d<p d(p' 

dr ~'~ dr' 

At the plane of symmetry 7i=a so that 
dr ldr\ , dr dr' 



= -i 



> ^^ rf7=; 



o?a: \dx I dy dy 

Therefore for r=r' we have 
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o da) dr 
dx dx 

That is, at the plane of symmetry the 
displacement perpendictilar to it is twice 
as great as that due to either wave acting 
separately, and the displacements paral- 
lel to this plane and the vertical displace- 
ments are equal to zero. Suppose now, 
further, that a^—o: then we have for 

d<p _ dcp ^ 

dr ~ dr' ' 
therefore, 

^dm dr ^ dq> 

dr dy dz 

From which we have — ^if there is no 
difference of phase between the waves 
from the parallel axes — at the plane of 
symmetry there is no displacement in 
the direction of the axis of X, L e, in the 
direction perpendicular to this plane; 
also that the displacement parallel to the 
plane and the vertical displacement are 
twice as great as they would be if there 
was but one wave. 

The reader who is interested in the 
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subject of wave motion will do well to 
read an article on the subject by Lord 
Bayleigh in the April number of the 
JPkilosophicai Magazine for 1876. An 
article in the September number of the 
same publication for 1878, though not 
bearing directly upon the subject, will 
also be found to contain much that is 
of value and interest; the article referred 
to is entitled " Hydrodynamic Problems 
in reference to the Theory of Ocean 
Currents," by M. Zoppritz. The mathe- 
matical theory of wave motion remains 
pretty much as Airy left it when he com- 
pleted his work on the subject — so no 
better reference can be given than to 
that work — ^for any one wishing to ac- 
quire a thorough knowledge of the sub- 
ject. 

§ 5. 

PBEE VORTEX MOTION. 

We have already seen under what cir- 
cumstances it is impossible for rotational 
motion to exist in a fluid mass. If the 
fluid in its initial condition has irrota- 
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tional motion— or, if it be at rest, and 
motion is induced by a system of con- 
servative forces, then the motion will 
always be irrotational ; i. e., if the quan- 
tity 

tutx+vdy-^iodz 

is at any time an exact differential it will 
always be one. The conditions for this 
quantity being an exact differential are 

dw do du dw „ 

dy dz ^ dz dx 

Suppose that these quantities are not 
equal to zero, but that we have 

P^l/dw dv\ 
__1 /du dw\ 

^"2 \di "■ dxr 

1 /dv du\ 

The quantities Sy 7 and 5> as is well 
known, denote the components of angu- 
lar velocity around the axes of jp, y and 
Zj respectively, of a particle whose ve- 
locities parallel to these axes are u, v, w. 
That these quantities should have certain 
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definite values different from zero is, of 
course, the conditioix that vortex, or rota- 
tional motion exist in the liquid. These 
values of S, t} and 5? pre-suppose that 
we know the values of </, v and to. A 
problem that now immediately presents 
itself for solution is to find the values of 
w, 13 and to, supposing ^, 7 and 5 to be 
given. 

Assume three functions tJ, V and W 
such that 

^ dy " dz 

_dV dW 
~" dz "" dx 

_dY dV 

^ dx^ dy^ 

The quantities t/, v and w must satisfy 
the equation 

du dv dw 
dx"^ dy ^ dz ^^' 

It is found without difficulty that this 
equation will only be satisfied by the 
above values of t/, v and to if the follow- 
ing conditions hold, 
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da^ + dy' + <?«*" ''' 
af'W <f'W (?'W _ „^ 

dx dy dz ~ ' 

The integrals of the first three of these 
equations are well knovm to be given by 

— /*/*/ S'dxfdy'dz' 

2^ J J J ^(iB'^) ' + (y '-y)* + (z'-zf ' 

V= 

1 /• /• /• tj'dx'dy'dz' 

^J 'f J V(^-*)' + (y'-y)* + (2-2)" 

W= 

'^ r r f Z'dx'dy'dz ' , 

where x','y,z' are the co-ordinates of any 
other point in the vortex element, and 
S,' fji' C are the angular velocities at this 
point; the denominator, 

r= V(^^-iY + {y'-yf + (i^^W~ 
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denotes the distance between this point 
and the assumed point to which the U,y, 
W refer. Before going farther it will be 
convenient to give two of Helmholtz's 
definitions. The line passing through 
any point and coinciding at all times in 
direction with the instantaneous axis of 
rotation of that point is called a vortex 
line. If we consider a number of vortex 
lines passing through every point in the 
perimeter of an infinitely small surface, 
they will cut from the rest of the fluid a 
filament which is called a vortex Jilament, 
or a vortex filament is an infinitely small 
filament of the fluid whose bounding sur- 
face is made up of vortex lines. 

Now, in our equations giving TJ, V, W, 
the points x, y, z and x' y', z' are sup- 
posed to lie on the same vortex filament ; 
we can represent an element of this fila- 
ment by dr, then our equations become, 



^= if'-?' 



tfdr 

r 
Vdr 

27t 
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where the integrations are of course ex- 
tended over all the space which is sup- 
posed to be filled with vortex filaments. 
Now to examine the condition 

dx "^ dy dz ■~^' 
By differentiation we have 
dJ} __!_//• f {x-x)5'dx'dy'dz ' 



_ _1^ /• /• /* {y-y')ri'dx'dy'dz' 
dy^'^'lnJ J J r* 

dW___J^ /• /• /' {zr-z')i:'dx'dy'dz' 

Integrating by parts we have, 
^-_1./ f S'dy'dz ' 

'dy" ^TtJ J r 
dz ~~2n<^ J r 
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Now since 

dS dtf dZ 

d^'^dp'^Jz^^ 
throughout the entire mass of the fluid 
we have 

dXJ dj dW^_ _1 j /• f B^'dy'dz' 
dx "^ dy ^ dz ~'^27t{*^ J \r 

This can readily be changed into a sur- 
face integral. If da denote an element 
of the surface of the vortex filament and 
cos a, cos /?, cos Yt denote the direction 
cosines of the outward normal to this 
surface, we have 

dx^d\/*=da cos y^ dx'dz'=d<s cob /?, 
dy'dz'^da Q,oB(i\ 
therefore our integral becomes 

taken over the entire surface. Now from 
the equation 



d^ dfj dz 
dx dy dz ""^ 



we have also 



ffA^+Ty + j^^^^y^^=o 
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by integration this becomes 
ffSdydz^ffrfdxdZ'^ffZdxdy=io 

or as a stirface integral, 

f{^S cos a + 77 cos /? + C cos y)d6=o 
consequently 

dx "^ dy"^ dz '"^* 

Substituting now the obtained values of 

<ZW dN 
TJ,V,W, in the equations ^=-> ""^ <^^- • 

and we obtain for w, v, to, the following 
values: 

-:'(*-;«:')]<?«'<?y'<?z' 

or, as these may be expressed, 
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Bepresentmg each of these differential 
expressions by w', «', w' respectively we 
see that u\ v\ v) are the increments of 
u, Vy w, which correspond to the element 
dx dy dz of the vortex fULament Writ- 
ing for ^convenience of reference the 
equations 

, <? /^* d^\ 

we see that they give rise to the equation 

This shows that, considering w', i?', w' as 
the components of a certain new velocity, 
the direction of the resultant 



of these components is at right angles 
angles to the direction of the axis of ro- 
tation of the element dr. Again, we 
have 
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,dr , dr , dr 
dx dy dz 

and the direction of this restdtant is also 
at right angles to the line r joining the 
element dx to any other. Thus we see 
that each rotating element of the fluid 
mass implies in every other element a 
Telocity whose direction is at right 
angles at the same time to the axis of ro- 
tation of the first and to the line joining 
the two elements— i.6., at right angles to 
the plane containing the second element 
and the axis of rotation of the first. It 
is easily shown that 

, dx ^ sin S^ 

when ^ denotes the angle between r and 
the axis of rotation. From this equation 
we see that the magnitude of this induced 
^velocity is directly proportional to the 
volume of the first element, its angular 
velocity and the sine of the angle between 
the line joining the two elements and 
the axis of rotation; and also inversely 
proportional to the square of the distance 
between the elements. Denote the an- 
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gular velocities at the time ^=o by S^, rf^y 
Co then the last equations of chapter I 
become 

^ j^dx dx . dx 
^^^Qa^'f^db-^^^dc 

"^^^'da^^'db^^'dc 

f,__pdz dz dz 
^"^'da'^'^'dh'^^'dc 

a, d, c being the co-ordinates of an 
arbitrary particle we have for the co- 
ordinates of another situated indefinitely 
near this a-\-da^ h-^-dh^ c-\-dc, and at 
the time t the co-ordinates will be x, y, z, 
x-^dx, y-\-dy, z+dz; now suppose that 
at t=o we have 

da_^dd_^dc 

that is da^ dh^ cd proportional to the 
initial angular velocities — and suppose^ 
further that the direction of this indefi- 
nitely small line coincides with that of 
the axis of rotation. Call the common 
value of these ratios £, € being an indefi- 
nitely small quantity, independent of the 
time ; then we have 
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da=S^€, db=zr^^e, dc:=z^B, 

Substituting these vaJues in the above 
equations, and we have 

da=S€y dy^TfB, dz^^S^ 
or 

dx__dj/^dz 

consequently the direction of the line 
joining the indefinitely near elements 
will at all times coincide with the direc- 
tion of the axis of rotation. This com- 
bined with our definition of a vortex line 
shows us that every particle of fluid that 
lies on a vortex line at any instant will 
always remain there. If we call go the 
resultant angular velocity we have 

^= a/5" + r^ + :=€Vdx* + dt/* + dz^ 

or the angular velocity so varies as to re- 
main always proportional to the distance 
between the two particles. We have all 
along supposed the density of the fluid 
equal to unity. Remembering now our 
definition of a vortex filament, we see 
that any vortex filament must remain 
composed of the same fluid particles. 
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CaUing now k the cross section of any 
filament and I an indefinitely small length 
of the filament, that is 

we have A:^= const.; but I is proportional 
to fl>, therefore fl>^= const, or the product 
of the angular velocity of an infinitely 
«mall portion of any filament into its 
cross section is a constant. Call now h^^ 
k^ the cross sections of a filament at points 
whose angular velocities are given by Wj, 
it>^ ; and let dr denote an element of the 
filament ; 

/, idB, drj d%\ 

=:^y d(T[^S cos (nx) + jf cos(ny) 

+ C cos (w«)] = — y dtfio cos (am). 

Now from the values of ^, 77, C, we see 
that the factor of dr in the left hand 
integral is=o; 

\fdGui cos (a>/l)=0 

But at the ends of the portion of the 
vortex filament that we are considering, 
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we have cos (o*n)=-i- 1, and for all other 
points cos {wn) ==a, and consequently our 
integral is equivalent to 

or the product of the angular velocity 
and the cross section is a constant 
throughout the vortex filament. As 
each rotating element of the fluids implies 
rotation in every other element, we have 
that all the particles of fluid must be in 
motion, and consequently from the defin- 
ition of vortex lines we see that these 
lines and consequently the vortex fila- 
ments cannot terminate within the fluids 
but must either terminate in its surface 
or must return into themselves; the 
former of these cases is illustrated by 
the vortices formed in running water,, 
and the latter by smoke rings. 
If in the expression 

J J J \ dx "*■ dy 

dz J 

we substitute for w, v and w their values 
in terms of the derivatives of U,V,W, we 
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/'Hf^-a-'-—' 



find that this is equal to zero, and from 
it we obtain the striking equation 

\dz '^ dx / 
jdY dV\ \ 
\dx " dyl J 

Idv du\ \ 



or 



But we may write for U, V, W the values 

But the expression for the energy of the 
fluid is 
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and by virtue of the above we have 

"We can now take up the simplest case of 
vortex motion viz., that in which the mo- 
tion is parallel to one plane. . If we as- 
sume this plane as xy and further make 
the motion independent of z we have the 
angular velocities around x, and y, and 
the velocity in direction of z equal to 
^ero, or 

We have thus 

~~ dy^ "*" dx^ ~~dx dy 
from which 

'■" dx^ ■*■ dy* 
which gives 

W=~/ciog/>«r(r 

when da is an element of the plane xy. 
Of course as C is independent of z we 
might have obtained this by integration 
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from our gcnenl faloe of W beiove 
^^cn. Vbtb p icpreBeBtB the <1iHtiiir8 
of the dement d4f in the plane xy from 
taxj cffhar point in tlist plane: Eedi 
Tortex filament implies in any otiier par- 
ticle of the fluid a Telocity whoee conn- 
ponentsare, 

1 Zdtf dp - 1 Zdff dp 

— . -TT- and — - . -r- 

n p ay n p dx 

and whose magnitade is 

1 ^da 
n p 

The direction of this Telocity is giyen by 

the cosines 

dp dp 

d%f dx 
and of the line p by 

dp dp 
dol dy 

or the direction of the velocity is at right 
angles to p. Assome two quantities a;., 
y,, which define by the equations 

x^f^da=zfxZd(r, 
^ yJZda=^fyZda. 

Now if we regard C as the density of a 
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mass distributed over the element dff of 
the plane xy we see that x^, y^ will rep- 
resent the co-ordinates of the center of 
gravity of this mass. Now /, the length 
of an indefinitely small portion of onr 
yortex*filament cannot alter, consequently 
by virtue of the equation A;^=: constant^ 
k or da cannot vary with the time and 
by virtue of A;a>= constant, <w or C cannot 
vary with respect to the time; conse- 
quently, we have, by differentiating with 
reference to t 

Substituting 

dt n^ p ' o' 

dy 1 /•:'( 

we have 
dx. 



P 
'da* x—x^ 
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The double integrals are=o therefore 

or the center of gravity of the filament 
does not change with the time. In the 
case of only one vortex filament let us 
write 

for particles as at finite distance from 
the filament we have 

dW dW ^^, 1 , 

for particles indefinitely near the fila- 
ment we see that W, n, v, are infinite and 
depend upon the cross section of the fila- 
ment and the angular velocity p. We 
also know that at the center of gravity of 
the filament u and v=o. Each particle of 
fluid that is at a finite distance from the 
filament we see has a uniform velocity of 

— and moves in a circle whose center is 

Ttp 

the center of gravity of the vortex fila- 
ment. Suppose we now assume a number 
of filaments whose cross section is indefi- 
nitely small. Write in general 
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and let Xi yi denote the co-ordinates of 
the centers of gravity of the filaments at 
the time t and pi their distances from 
the point (icy). Then for all points at 
finite distances from the filaments we 
have as before 

Now each filament inducing a certain 
amount of motion in every other particle 
of the fluid induces motion in the centers 
of gravity of every other filament, there- 
fore the filaments change their places in 
the fluid. But here, as before, that por- 
tion of w, V which each vortex filament 
gives to its own center of gravity is =o. 
Suppose that the point from which the 
ps are measured is one of the centers of 
gravity, for example x y. This will ma- 
terially simplify the investigation by con- 
fining us exclusively to the influence of 
the system of vortices upon its diflPerent 
members and as this point a?, y, is an ar- 
bitrary point no generality is lost. We 
have thus 
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dW, rfW, 

"• dy. ''•- dx, 

(m, log /)„+m,log />„+ . . .milogpli) 

or briefly 

1 i 
W,= — - -2 Wilog/t>w. 

We can now assume a function Q such 
that 

Q= — -2m m log p 

and then we have 

dx. dQ .cTaj, dQ 
"^^ dt "dy^'^^dt-dy: 

"^^ dt ^"d^^'^'dt^^d^: " " 

A complete system of integrals cannot be 
in general obtained, but by observing one 
peculiariiy of Q we can obtain two inte- 
grals — 

P ij=V{Xi — xjY-^-ij/i — yj y 

If we increase Xi^xj^ or yi,yj,hj the same 
quantity, jOi^will be unaltered; this gives 
lis then 
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dxi ' dt/i 

or d^i dxi 

i at i at 

from which 

2mi Xi = const. 

2mi yt = const. 

or the center of gravity of the system of 
Yortex filaments is unaltered. Again 
since 

we have 

c?Q=o, .-. Q= const. 

the equation of the lines of flow of the 
fluid. Introduce now polar co-ordinates> 

x^—p^ cos 3", x^=p^ cos 3^, -h ... . 
y,=/>, sin 5, y,=/?, sin S", 4- 

we have by these substitutions 

dp^ __ dQ, dp^ __ dQ 

^'^' dt "■ d§f: "^^^^ di " d^ 



d^^ dQ dSr dQ 

""'^^^-dt =''dp:'"^^^-^t ^"Tp, •• 

If now we increase all the ^'s by the same 
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quantity, Q will evidently remain unal- 
tered and we have the equation, 

The first row now gives by addition 

^^ ^nii Pi *= const. 

Now let us suppose d to remain un- 
changed but p to become np, then log p 
becomes log p + log n, and log pij be- 
comes log Pij 1+log n, and in conse- 
quence Q will be increased by 

"" "Z'L^i^, log ^ + m^m^ log n + m^m^ logn 

+ . . . ,mj Tij log w+] 
or 

log n 2 mi nij 

and consequently we have 

^ ^Q 1 ^ 



^ log Pi n 

or 

^pi "f— = 2 nii mi 

clpi n - ' 
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But we have also 



_=«^,^,„ 
Substituting this value gives 

2 mi pid^i = — Snii nij 

7t 

Assume now the case of only these 
vortex filaments existing in the fluid. 

The equations ZtrnXi = const. ; and 
2mii/i = const, become 

m jPjC0s5j + w,p,cos5, + m,/o,cos5, = C, 
m^pfimB^^ + m,p,sinS^, + m^p^&mS^^=C^ 

Multiply the fir^t equation by cos S^, the 
second by sin S^ and add, 

m,/Oj +m,/o,cos(S',— S^,) + m,/o,cos(S^,-j5) 

=C,cosS', + C,sinS', 

Again multiply the first by sin ^9 and the 
second by cos S^, and add 

»7i^,sin(S^,— 5,) +m^,sin(S^,— S^J 

=0,cosp3—C, sins', 
Again, 

Q = -~ [^i^Jog/o,, + w,w,logp,, 

+ 7w,m^log^„+ ]=const. 

and 

m^p* + WaP,' + ^3/}/= const 
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Through these four equations we may 
express any four of the quantities p^^p^tPa 
S'.-S',, ^^—^^ in terms of the fifth ; for 
example p, ; then the equations 

" ' ' dp^ 
and 

3 J^ 3 

1 7t I 

will enable us to express S^^ and t as func- 
tions of p, and afford a complete solu- 
tion of the problem. 

Assiune now only two vortices, and 
take the origin of co-ordinates at their 
common center of gravity. This point 
do5s not move, and we have 

'dt " dt 
Q in this case becomes 

——m^m^ log (p, + pj =con8t. 

and also 

2miPi = niiPi 4- niiPi = const, 
from these two equations we obtain 

Pj = const., p, =const. 
Again the equation 
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becomes 

/// 

giving 

d^^ __ dS^^ ^ 1 m^m^ 

'dt ~~dt '^ 7t' m^p*-\-mj)* 

If the direction of rotation of both 
vortex filaments is the same m^m, [which 
depend on 5, and 5 J have the same sign. 
But suppose m^=^m^, then 

dS^^ 1 m, 
dt "" 7tp^'-p,^ 

But we have now for the center of 
gravity 

y__ — i^j »^«-=oo 

or the center of gravity of the two fila- 
ments lies at infinity. Their velocities = 
their angular velocity x by the distance 
from the center of gravity differ from 
each other by an infinitely small quantity 
and can be expressed by 
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2 dt 

but by our preceding equation giving the 

value of -TT this is 
at 

1 m. 



2^ />,-/>, 
the direction of the motion is, of course, 
perpendicular to the line giving the cen- 
ters of gravity of the two filaments. The 
particles of fluid lying between the fila- 
ments move forwards in the same direc- 
tion as do the filaments, the one-half way 
between them moving four times as fast. 
Let us suppose that the vortex filaments 
at the beginning of the motion lie on the 
axis of X at equal distances from the 
origin, then the particle above referred 
to will lie at the origin. Also write 

^ — —=a the absolute distance of each 
z 

filament from the origin. We have then 

for the co-ordinates of the filaments at 

the time t {a, y') and, by virtue of what 

has been said the co-ordinates of the 

particle originally at the origin will be 
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(0,4y'). The equations of the lines join- 
ing these two points are 

«> y, 1 

o, 4y', 1 =0 

±a, y\ 1 

The intersections of these lines with y=o 
are given by 

«= ±y or «= ±f (p^-p^) 

that is, the lines joining the particle half 
way between the filaments with the cen- 
ters of gravity of the same pass through 
fixed points on the line joining the origi- 
nal positions of the centers of gravity of 
the filaments. The points he outside of 
the original positions of these centers of 
gravity and at an absolute distance from 

them=J(p,-/o.)=|, ^9^'. 

Tit 

The particles of fluid that lie in the 
plane bisects at right angles the line join- 
ing the two vortex filaments will remain 
in this plane. If this plane be considered 
as a fixed boundary, we have by con- 
sidering one of the filaments the case of 
a filament moving parallel to a fixed 
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plane which limits the extent of the 
fluid. 

Let us now assume that the vortex 
filaments are so arranged as to form the 
continuous surface of an elliptic cylinder 
of finite cross section, and further 
assume S, as constant for every point of 
this cross section. As the same particles 
of fluid constantly remain in any vortex 
filament, the bounding ellipse of the 
cross section of the fltiid will always be 
composed of the same fluid particles. 
The equation of this line will be a func- 
tion of a;, y and U £^d may be written 
for the moment as 

/(a!,Y,0=a 

Then by virtue of the above we have 
generally 



or 



But, 



dt dx dt dy dt 



df ^ df ^ df 
dt dx dy 

dy^ dx' 
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Therefore this equation becomes 

df dW^^dW df_^ 
dt dy dx dx dy 

Now the general equation of our ellipse 
is 

when a, )^, y are functions of t Assume 
another system of co-ordinates coinciding 
with the axes of the ellipse and also pass- 
ing through its center; then 

a;'=a5Cos 0+y sin 6 
y'^z-^x sin O+y cos d. 

Call a and h the semi-axes of the ellipse 
then we have for its equation, 

a' ^ b' •"^• 
Substituting for «', y', their values as 
given in terms of x and y, and this 
becomes, 

, ^>'cos'S'-{-a *sin'S' , „ [{b* -«*)cosS'sinJ^i 

^ S^P ^^ ^^ ^ 

. , ^'sin'S^ + a'cos'S' ^ 

-"^ ^^ =1 

Comparing this with 



and we obtain, 

a^by^b^Bm^^-^-a^cos'^ 

In these a and b are constant, but a, /3y 
y^ ^ are functions of t Now W satisfies 
the equation 

for all points in the interior of the ellipse, 
and its value is obtained by integration 
of the equation 



00^ «» y' 






. d\ 



«K^iog:^,+i)-^,(^'*+-0 



The integral of this is 

2 

or, since we only use the derivatives of 
W, we may write it, 

for all interior points, and for points at 
the boundary. If now we write for 
brevity, 






145 

A.=h cos' S'+a sin* *& 

B=(ft— a) cos S^ sin S^ 

r=6sin'3^+acos*J^ 
we have 

Let ns examine again the condition that 
we obtained for the bounding line of the 
cross section of tije cylinder, yiz: 



dt dx dy dy dx ~ 



we have 



df da ^ ^d/3 d^ , 
dt^W-^^-dt"^-^ dt^ 



dy dx ^ "' ' a+b 

Equating to zero the co-efficients of a;' xy, 
and y* separately we have 

(a+*) ^=45(aB-/JA), 
(a +6) ^=45(/Jr-yB). 
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If ^ cannot be determined as such a 
function of t as to satisfy these equations 
then will our equation of condition 

df dfdWdfdW_ 
dt dx dy dy dx " 

hold for all points in the cross section. 
Forming the derivatives of a, /3, y with 
respect to t and transforming them by 
reference to the values of A, B, P we 
find that the function of t sought is given 
by the equation 

d^ ah 



dt {a + by 

We have thus the value of the angular 
velocity with which the cylinder rotates 
around its axis. The rotation of the 
cylinder also induces relative motions 
among the component vortex filaments. 
These are obtained by regarding x' and 
y' as f imctions of t We have by differen- 
tiation 

~dt "^ ^ dt dt^ dt 

and the other components of the veloci- 
ties in the directions of x' and y* are 
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dW dW 
dy'' "^ dx' 

Therefore we have by combining these 

dx^ _dW ,d^ 
dt "" df ^ dt 



dy[ ___^ _ , d^ 
dt "" dx' dt 



Now, 



and 

d^^ ^ ab 

dt"^"^ {a-\-bY 
Therefore 

dx'_ 2ga' , dy'_ 2gy , 
Differentiate each of these for t and with 

and we have after integration of the 
resulting well-known powers 

a;'=a;cos {dt+i) 
y'=b I cos (6t + i) 

when I and i are the constants of inte- 
gration and determine the particle of 
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fluid to which x* y' have reference. p> 
1 because then for the cases when 
cos (0^^-^)=l we should have a;'> a 
which cannot be, .*. ^ is a proper frac- 
tion, d of course denotes the angular 
velocity of the cylinder, or 

from which t?= Qu 

Solving the equations 

xz=.x cos 5 + y sin 5 
y'z^. —a sin S" -h y cos S" 

for X and y gives 

x-='X' cos S' — y' sin S^ 
y=a/ sin S^ + y' cos 5 

Substituting for »', y' and J& their values 
these become 

25= aZ cos (dt-Vi) cos 0^ 

-Wsin((9< + i) sin (9^ 

y=a? cos (0^^-^) sin dt 

+ 5/ sin (dt-\-i) cos 0^ 

by expanding the quantities cos sin 
{Bt'\-i) and collecting the terms these 
equations m9.y be written 
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x= -J^ I cos {Wt^i) +^ cos i 

y= ""l^l^m (2et-^i)-'^l sin i 

by differentiation with respect to t we 
obtain 

dt 



^=(a-^b)ei COB {2et + i) 
at 



ab 



dy 
dt 
from which 

and also, 

(x-'^ I cos i)\[y + '^ I «^i) = 

From these equations we see that each 
vortex filament moves in a circle with 
uniform angular velocity, the tii3[ie of ro- 

tation being evidently 3-; and that the 

position of the center and the radius of 
the circle varies for different filaments. 
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Suppose one of the semi-axes a or ^ to be 
infinitely greater than the other; this 
gives 6^0, and eonsequently S^=o, or 
the straight line which has become the 
limiting case of the ellipse does not 
rotate. If a==^ our ellipse becomes a 

circle, and we have ^= ^ 5 ^ *^s case 

there is no change of the relative posi- 
tions of the vortex filaments, but they 
all rotate around the common central 
axis with angular velocity 5- 

The next case of vortex motion that 
we shall consider is that in which the 
vortex lines are circles having their cen- 
ters in the axis of z. The direction of 
the axis of rotation of each fluid particle 
will lie in a plane at right angles to z and 
be parallel to the tangent to the. vortex 
line at that point.. The reader will do 
well to observe the motion of a ring of 
tobacco smoke; he will see that the ring 
seems to be turning inside out, each 
particle moving in a plane passing 
through the axis of the ring and revolv- 
ing in a circle whose axis is in the direc- 
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tion of the tangent to the ring. He thU 
also obsjaorve that there is no motion 
around the luxis of the ring. Now if- we 
introduce pol&r oo-ordiBates p and S^ we 
have, evidentij, 

xcz/j cos S' y==:/> einS^ 

and for the rotations we may obviously 
write 

5=— A sin S', t^=X cos 5, C=o, 

when A. is not a function of S^. The 
equation of the path of the particle is 
evidently one between a, y, z where x and 
y are connected by the relation 



so that the equation of the path can be 
made to depend only upon p and z* Re- 
suming now our equations 

_dW _ cIY _cnj ^cTW 
"^ dy dz^ ^ dz dx^ 



and observing that 

W=— / — 
27n^ r 

ttiese give 



_dY dV 
"^dx "" dy 
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_ <?V dV ^dY dV 

dz dz dx dy 

If now we assume an element dt* given 
by the coordinates S^', />', z' and for 
which \ becomes V\ also denote by r its 
distance from the point (p, S^, z) or (aj y g). 
Now onr equations for U and V are 

V sin ^dr' 






r 
V COS :&'e?r' 



2;r«^ /• 2;r«^ r 

When 

(?r'=efo5' dy d^'^p'dp' dz' d^' 
and 

^= ^(2-3)"+ p' + p'*-2pp' cos(5'-5) 
Now make for convenience 

this gives S^'=^+S'and^S^'=<?9> there- 
fore we have again for TJ and V, 

A^[sin y cos S^-f-cos^ sin ^'\dqi 
V{z'-z) • + p* + p'* — ipp' cos 9> 
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^= L/f^'^^'/^^f 



2;r« 

A' [sin q) sin S^— cos ( p cos ^ ']d(p 

Viz^zy + p' + p'^'-^pp' COB ^ 

In both cases we have the integral 
/» sin ^ dq) 

•^ \/(2'-s)»H- p*+ p'*--2p/>'cos 9? 

this is to be taken between o and 2;r, 
these being the limits of S^'— S', 

/* sin <p d q) 

2;r 



=log VC^'-^)* +73' + p''-2p/o'cos9?] 







— 1 (g^— g)* + /Q* 4- />^'- 2pp'coQ(p _ 

-2 ^^^ (g'-g)* + p'' + p'*-2/Q/9'cos9>~^ 

and our expressions are thus reduced to 
/• A'cos^sinS'c?^ 

•^ A/(g'-g)' + />' + p"-2/9/9'C0S^ 

/» A' COS 9? COS 5d^93 

*^ V'(g'-g)" + /o'+p"-2/}p'cos 9> 
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Denoting the common integral by ^ we 
have 

U= - ^/p'dp'/dz' ;i' # sin S 

V = ^fp'dp'fdz* V 9 cos Sf 
^ will clearly be of the form 

XJ and V now differ only by a factor, in 
fact we have 

V=-Utan^. 
So if we write 

we will have briefly 

U=-Psin:^, 
V=3 P COS ^. 

The value of the function ^ is not diffi- 
cult to obtain ; we have 

^= f cosy^y 

J^ V(2'— 2)' + p' + /9'*— 2pp'C0S9 

We will find the integration much sim- 
plified by the introduction of a new vari- 
able W defined by the equation 
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2 

Then we have 

cos 9?= -[1-2 sin V] 
cog ^<p=f2[l-2 sinV]€^i^ 
and for the limits we have when tp 
=(o,2;r),^=(^,-Z,| Making this 
transformation we have 

"2 

^=2 / ^ Cl-2 sinV1(^ 



1/ a/(^'-. 



-2)' + p' + /o'''+2/}p'-4p/osinV 
or, 



-1/? 



[l-28in*y]jy 



a/(s'-s)' + (p+p')'-4/o/o' cos 9 





V(3-a)'+ (p'+p)Vi . 4pp'sux S*- 

• ^r I    * - I 
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+4/? 

2 sin*^ d^ 



V(2'-z) » + (p' + p) » ^ _ 4p/3'Bm'y 
Make, 

then 

1 A: 



then 

2 

2* /• t7^ 





n 

2 

2A? /•_2sin^M^ 



The first of these is the complete elliptic 
integral of the first kind; we will, as 
usual, denote it by K. Examine now the 
second integral ; we hare on multiplying 
it numerator and denominator by A; 
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n 

2 

4 p Tc^ Bva}(pd^ 



-,/ 



^^fpp'*^ \/\^k^ sin^ W 



n 

2 

,l-(l,A?»sinV) 



* pX.—\X.—K bill 9^^ 



71 

2 
4 /• 6?^ 



./ 



^^/pPo Vl-/<j' sin ' ^ 



n 

2 
4 



The second of these is the elliptic inte- 
gral of the second kind and denoted by 
E; we have then finally for ^, 

_ — K 4 4 

or 
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and conseqnently 

sin ^ dp'dz\ 

COS ^ dpldz\ 

In thjB function that we have denoted by 
^ we see that the derivatiyes taken for 
z and z' have the same absolute values 
but opposite signs; consequently 

when da=dpdz. But we have also 

dV 1^ y*-^, ^^ , , 

clz'^ 2;r«^ ^ dz 
therefore 

Jp-XdiS^o, 

But we have for the energy T the equa- 
tion 

T=/^7r(U5+V7) 
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Substituting for TJiV^S and 7 their values 
this becomes, 

integrating with respect to ^ from to 
2 ft, and writing again dpck^diX, 

T=27rf'PpXd<T. 

Substituting for P its value we have 

Tl=if$fip'XVdcfd<f' 

when of course dcf and diS' denote the 
cross sections of the vortex filaments 
under consideration. Let S denote the 
component of velocity in which p in- 
creases : 

S=Vw*+v* 

which is the same a^ • 

w=S cos ^ t;=S sin ^, 

But 

dN . dV ^ 

dz dz 

y= — = sm ^ 

dz dz 



therefore 

8=-: 

dz 



s=-li'. 



Now for w we have 
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dx dy 

and substituting for U and V their 
values. 
This gives 

dp p 

dp s 

and we may also write 

*/,=-^) since ^=0. 
dz dz 

yFrom these we have the equation 
^ dz ^^~^ ^ *^® iorm 

J p8\d<5^=.0 

and also 



therefore 



dt dt 



fp'^^Xda^o 



dt 

and since for each vortex filament Xdff is 
constant this gives 

//0*Ada'=const. 



161 

Some other interesting forms may be 
given before we proceed to the examina- 
tion of a speciaJ case. We had 



k^ = 



^P'P 



Taking logarithms this becomes. 

2 log ^=log pp'-log K^'-a:)* + (p' + p)>] 

Differentiating with respect to p and z we 
obtain 

k^dp''{z'-z)^ + {p'-\-pY 
2 elk zz(z'-z) 



k dz (2'--2)*-h(/o'+/9)* 
consequently, 

kYdp ^dz\ (2'-2)*-|-(p'+p)* 

The denominator of this second member 
does not change by the interchange of 
accented and imaccented letters, but the 
numerator does change its sign, also k 
does not change by making the same 
transfer, therefore 

dk t^dk 
dp dz 
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assumes the opposite value by writing 
the accented letters in the place of the 
unaccented and vice versa. If from the 
value of ^ in terms of K and E we obtain 

-jT we will see that this quantity does not 

alter by the interchange of accented and 

unaccented letters, consequently the 

quantity 

d^l dk dk\ 

dk\ dp dzJ 

assumes the opposite value after the in- 
terchange. We have by partial differen 
tiation of ^ 

d^^^d^ d^dk 
dp "^ dp dkdp 

_ c7a> ^ ^ 1 £ 
~ dk' dp'^ z p 

d0_^ d0 dk 
dz ~" dk dz 

Consequently the above quantity assumes 
the form 

d^ ' d$ 1^ 

dp dz z 

And by virtue of the property proved for 
this quantity we know that 
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But we have 
therefore 

Since now 

^ dp~ dp 

and ^_ 

(f(P/j) , <?I* 

this equation becomes 

f{wp-zs)pKda- - /vpXda=o 

T 

^^ f{wp-^s)pXda=^—^ 

We will now introduce the comple 
mentary modulus A;' defined by 

we have for this modulus 
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If k is very nearly equal to unity, that is, 

if h! is indefinitely small we see that ^ 

"will be of the same order of magnitude 

as K, and again, that P will be of the 

same order as #. We will examine the 

order of K on the supposition that T<f is 

infinitely small. Since Jd is indefinitely 

small, we have, neglecting higher powers 
than the second, 

z 
Now 

2 



^=/ 



a/I— ^*sin*^ 



=/ 



7t. 

2" 

d4> 





7t 



Vl— sin*v^ + A;'* sin*^ 



=/ 



dip 



V cos* V' 4- ^'^ sin" 

or, by introducing an indefinitely small 
quantity € which is, nevertheless, indefi- 
nitely large as regards A;', 
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n 
2 






^ ____ 



n 



f- 



V^os^^+A;* sin *^ 

Write now in the first integral ^ - in 

in place of W\ since throughout the in- 
tegral Q is small the integral becomes 

dd 



/ 



v^'^+k^o 

=A^^^ ^^^ — - 

or since A?' is infinitely small with regard 
to ke this is=| log ^ or = log^' 

In the second integral k' sin <p is 
throughout smaJl as regards cos <p and 
this integral is 
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2 



5-£ 



/— ^= log tan ] o^^-o^ f 
0*^ COS ^ (2 2 J 

or what is the same thing = log -. 

Hence we have 

2£ 2 4 

K=log-^ + log- = log ^, 

Consequently for k' indefinitely small we 
see that K is indefinitely large, and ^ 
and therefore P are indefinitely large of 
the order log k\ Of course ^ does not de- 
pend on E as for k nearly equal to unity 
we have, 

27t 



E= y cos (p d4»^l. 



Representing the elliptic integral of the 
first kind as is usual by F, we have 

<P 

d(p 



"=/; 



Vl— Ar*sin*^ 
and also 

2 

E=y(l-^asinV) #. 



167 

Differentiating for k^ and writing 
Vl— A?*Bin *^= J 

c?F_ pJc sin*^c?^ 

c?E_._ /•/bsin^^c?^ 

Now writing 8in*^=: -^ (1—^*) we see 
that these two integrals depend on 

the two first of these are F and E 
respectively; as regards the third, we 
have 

c? sin ^ cos ^__1— 2sin'^H-A;*sin*^ 

or 

,2 cZ sin ^ cos^ _J^— A;'^_^_ k'^ 

dip' A It 2"' 

and thence by integration, 

/^V^— 1 -pj^K^sin^cos^ 

The expressions for -^ , and -^r-, thus 
become, 
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and for the complete functions when 



dk 

dE 
dk 



^ 2' 

When k' is indefinitely small the first of 
these is of the order jj^^ and the second 
by k\ Now we have obyionsly if A?' be 
indefinitely small that-^ is of the order 
,,2. We had also 



dk_k (e'-«)»4./»-^a 



dp ^piz'-zY^-ip'+pY 
These quantities are of the same order 
as k^. Therefore 

d^ d£dkl ^ 
dp dk dp z p 
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and cW_Ji^ dk 

dz ^dk dz 

are of the same order as jy. By the aid 

of these preliminary investigations we 
will now proceed to the examination of 
the case when only one vortex ring exists 
in the fluid, and will further more sup- 
pose this ring to be of indefinitely small 
cross section and of the same order of 
magnitude as the indefinitely small 
quantity «. We may again write as be- 
fore 

m^=xf\d(X 

as m will be finite and as d(J is of the 

order f, A must be of the order -5- As- 

sume the equation of a circle such that 
the fluid elements of which it is composed 
shall lie indefinitely near the vortex fila- 
ment. Let its equations be 

We had 
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for all points lying at a finite distance 
from the circle. For these points by aid 
of the equation^ 

dz dp 

we can find the values of 8 and w. But a 
difficulty exists inasmuch as p^ and z^ are 
functions of the time and as such will 
have to bo determined^ and to that end> 
it is necessary to consider the points that 
lie indefinitely near the circle, or points 
in the vortex filament. 

Suppose that the two circles (/>, z), (j>^, z^) 
are at a distance apatt, that is of the same 
order as e. If we call r the distance 
between them, we have 

r8=/o«4-iO,«-2/>/)„cosB'+(2,-2j)« 
or since 3^ is indefinitely small 

this is of the same order as k^^ therefore 
k' is of the same order as a — ^therefore 
by our preceding investigations we see 
that P is of the same order as $ when k' 
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is of the same order as e. And also that 
P is of the order log. e. We had for the 
energy 

therefore T is of the same order as log. €, 
The preceding investigations taken in 
connection with the equations 

dz dp 

show that inside the vortex ring s and w 

are of the same order as-** We will now 

€ 

examine more closely p^ and 2,. We have 

zJp^^^d(T=fzp^\d(T 

We assume that A is constancy of the 
same sign, consequently the circle (/0^,2J«) 
lies either in or indefinitely near the 
vortex filament. 

Now we found that, \d(T not varying 
with the time, 

y /o*Ae?a'=const. 

consequently p^ does not vary with the 
time. That is, if only one vortex fila- 
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ment exists in the fltiid its radius will re- 
main unaltered, during tlie motion. 
Now to examine ^^, We have the 
equation 

JXd(S=^m 
from these 

Differentiating with respect to t 

But we had also 
^T= /(wp-ez)pXdff 

consequently 

^" eft ^n ^ ^ dx 

T, we have seen, is constant and infi- 
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nitely great, of the same order as log e; 
the difference between the values of 2, in 
the second member, we know to be of 
the same order as e, and we have seen 

that -/- is of the order of — 5 *conse- 
at e 

quently the second member of the 

right hand side of the equation is finite. 

Therefore, it follows that -^* is infinitely 

great of the order oi log e, and, neglect- 
ing the finite term — is constant. Since 

dz 
T is positive, -=^ must have the same 

sign as m, ^. e., as A. Thus, if only one 
vortex ring exist in the fluid, it will re- 
tain its radius unaltered during the mo- 
tion, and will advance in the direction of 

dz 
the axis of z lyith the velocity -^*. Now, 

as this vortex motion implies motion in 
all the particles of the fluid we have, that 
all the fluid particles at a finite distance 
from the filament flow through the ring 
in the direction of z, or the reverse, ac- 
cording to the sign A. If S^=o, we have 
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A =77, and, according to the con- 
vention for positive rotation we have, 
that the motion in the direction of z will 
be positive, if, in the case 0=0, 7 is pos- 
itive. • Therefore, it follows that the ring 
is moving in the same direction as the 
fluid particles are flowing. I will now 
give the concluding remarks of Helm- 
hotz's great memoir as nearly as may be. 
We can now readily see in general how 
two vortex rings having the same axis 
will move with reference to each other, 
by observing that each will have its mo- 
tion modified, due by the motion of the 
particles of fluid, due to the rotation of 
the other. Suppose that both rings 
have the same direction of rotation, then 
they will both more forward in the same 
direction, and the former will widen and 
move more slowly, while the latter will 
contract and move forward more rapidly, 
finally overtaking ajid passing through 
the former, when the same operation will 
be repeated, the rings continually chang- 
ing position throughout the motion. 
Suppose that the vortex rings have 
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equal radii, the result is not changed in 
the case of the same direction of rota- 
tion existing for both. But now let 
them equal radii, and equal but opposite 
angular velocities; they will approach 
each other, and both will expand — ap- 
proaching very near the effect of one 
upon the other is greatly increased — 
and they expand with constantly increas* 
ing velocity. 

Now suppose that the rings having 
equal .and opposite angular velocities, 
are symmetrical to each other. Then 
the motion in the direction of the 
axis of those particles that lie midway 
between the rings, is (?. We can con-r 
ceive this surface in which these particles 
He as a fixed boundary, and we have the 
case of vortex rings moving in contact 
with a fixed surface. These rings can be 
readily formed in water; or, rather, half 
rings can be formed, if we draw through 
the water rapidly, and for a short dis- 
tance, a . half immersed hemispherical 
vessel. Half rings will be formed in the 
water, having their axes in the fine sur- 
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faoe of the fluid, which will more exactly 
as described in the theory. The free sur- 
face of the watei* will form a limiting 
plane, passing through the axes of. the 
rings, and will not affect the motion. 
Rings of tobacco smoke have a rapid mo- 
tion forwards in the. direction of and due 
ta the impulfiivie force which prodtrbed 
them;' at the sanie time the ring flows 
through itself in the direction of the mo- 
*tion of translaition. 

li is very interesting to observe the 
motions of smoke rings, and for this 
pui^ose the foHowing simple apparatus, 
which has been described in a great many 
pltices, Will %e found useful: A rough 
box, abotit tei!» inches long, and the same 
height and width, is large enough; one 
end of the box to be open, and over this 
stitch a piece of cloth or rubber'; make a 
hole, about three inches in diiameter, in 
the opposite end of the "box, and a num- 
ber of slidefe having smaller holes in 
them, to be placfed ovei* th^ larger* open- 
ing and concentric with it. NdW placje 
inside of the box a vessel containing salt, 
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on which pour strong sulphuric acidi 
and also place in the box a piece of cot- 
ton saturated with ammonia; fumes of 
ammonium chloride will immediately fill 
the box. Now tap on the stretched mem- 
brane; rings will issue from the hole in 
the slide at the opposite end, and will 
move forward with velocities proportional 
to the force of the blow struck. A very 
light tap is all that is p.eGes8ary, and, in- 
deed, is all that can be given, if it is de- 
sired to investigate.the motion, as, other- 
wise, the rings move forward with such 
velocity that they can scarcely be fol- 
lowed with the eya If the rings are al- 
lowed to impinge upon a surface, the ro- 
tational velocity i^ suddenly increased 
very much, a;nd the rings thus spread out 
over the surface* 

The same effects will be noticed if two 
rings be allowed to meet each other in 
their motion through the air. If the ori- 
fice be elliptic, the rings will be seen to 
interchange rapidly their axes, vibrating 
about a mean circular position. 

If bubbles of phosphuretted hydrogen 
be aUowed to escape into the air, each 



